NEW  YORK  uMi:/L:;:.iTY 

INSTTTUTEOFMA 

L-. 

g»  Ml>^ht  ^i««r  *W   •'«»<*'  hH  t 


AEC  Computing  and 
Applied  Mathematics  Center 


AEC  RESEARCH  AND  DEVELOPMENT  REPORT 


PHYSICS 
TID-4500 
I4th  Ed. 


\ 


NYO-2877 


EXPERIMENTS  ON  SUPERSONIC  PLASMA 
FLOW  ALONG  MAGNETIC  FIELDS.  II. 

by 

David  M,  Wetstone 
September  30,  1959 


Institute  of  Mathematical  Sciences 


NEW   YORK    UNIVERSITY 

NEW    YORK,    NEW    YORK 


This  report  was  prepared  as  an  account  of  Government  sponsored  work.  Neither 
the  United  States,  nor  the  Commission,  nor  any  person  acting  on  behalf  of  the 
Commission: 

A.  Makes  any  warranty  or  representation,  express  or  implied,  with  respect  to 
the  accuracy,  completeness,  or  usefulness  of  the  information  contained  in 
this  report,  or  that  the  use  of  any  information,  apparatus,  method,  or 
process  disclosed  in  this  report  may  not  infringe  privately  owned  rights;  or 

B.  Assumes  any  liabilities  with  respect  to  the  use  of,  or  for  damages  resulting 
from  the- use  of  any  information,  apparatus,  method,  or  process  disclosed 
in  this  report. 


As  used  in  the  above,  "person  acting  on  behalf  of  the  Commission"  includes 
flny  employee  or  contractor  of  the  Commission,  or  employee  of  such  contractor, 
/o  the  extent  that  such  employee  or  contractor  of  the  Commission,  or  employee 
Jof  such  contractor  prepares,  disseminates,  or  provides  access  to,  any  information 

I  pursuant  to  his  employment  or  contract  with  the  Commission,  or  his  employment 

Iwith^uch  contractor. 


UNCLASSIFIED 

AEC  COMPUTING  AND  APPLIED  MATHEMATICS  CENTER 
Institute  of  Mathematical  Sciences 
New  York  University 


PHYSICS  NYO-2877 

TID-4500 

l4th  Ed. 

EXPERIMENTS  ON  SUPERSONIC  PLASMA 
FLOW  ALONG  MAGNETIC  FIELDS.  II. 

by 
David  M.  Wetstone 
September  30,  1959 


-  1  - 
UNCLASSIFIED 


NY0=2877 
PHYSICS 


ABSTRACT 

Time=lntegrated  photographs  are  presented  of  col- 
limated  plasma  flow  parallel'  to  magnetic  fields  up  to 
5  weber/m  ,  into  and  out  of  magnetic  mirrors,  and  along 
curved  fields.,   The  material  projected  along  the  field 
lines  was  substantially  Cu  plasma,  created  by  short 
vacuum  spark  bui*sts  (plasmoids) ,  and  injected  into  a 
vacuum  chamber.   Non-uniform  solenoids  expanded  or  com- 
pressed the  plasma  flow,  which  followed  the  field  lines. 
In  the  case  of  curved  fields,  more  complicated  flows  were 
observed,  due  to  the  presence  of  internal  polarization 
fields,  and  the  plasmoid  followed  the  field  lines  only 
when  the  density  (up  to  3  x  10^  ions/cc)  was  sufficiently 
reduced,   Photomultlplier  studies  were  employed  with  some 
geometries  to  estimate  center-of-mass  and  expansion 
velocities  which  reached  5.5  and  1.9  cm/jxsec,  respectively, 
adequately  demonstrating  supersonic  flow.   The  latter  figure 
gives  an  upper  bound  to  ion  temperature  of  120  eV  for  Cu, 
This  ion  energy  is  considerably  higher  than  those  normally 
encountered  in  spark  channels  and  is  a  substantial  fraction 
(0,12)  of  the  axial  energy,   A  mechanism  is  described  for 
achieving  such  thermal  energies  through  the  interaction  of 
a  strong  shock  wave  with  the  initial  discharge  channel. 


-  2 


NYO-2877 
PHYSICS 


thermalizing  a  large  fraction  of  the  axial  energy.  The 
upper  bound  for  ion  Larmor  radius  was  that  of  the  plasma 
flow.   Estimates  are  presented  of  electron  temperature  and 
degree  of  ionization.   A  theoretical  total  reflection 
coefficient  is  derived  for  a  plasmold  entering  a  magnetic 
mirror.   An  experimental  transverse  drift  velocity  of 
5  X  10  em/sec  (in  the  straight  solenoid)  is  observed 
and  discussed.   In  the  Appendices  a  complete  description 
is  included  of  the  experimental  equipment  and  electrical 
circuits,  and  a  digital  machine  computation  method  is  des- 
cribed for  obtaining  field  line  distribution  and  density 
in  a  wide  variety  of  coll  geometries. 


-  3  - 


NYO-2877 
PERSONNEL  AND  ACKNOWLEDGMENTS 

This  project,  a  part  of  the  Project  Sherwood  program  at 
New  York  University,  was  under  the  general  supervision  of 
Professor  Harold  Grad,  Institute  of  Mathematical  Sciences,  and 
extended  over  the  period  of  February  I956  through  June  1959. 

This  present,  and  final,  report  deals  largely  with  work 
performed  since  April  I958. 

The  project  was  conceived  and  initiated  by  Professor  David 
Plnkelstein,  who  was  its  director  until  July  1958^  his  contri- 
butions to  the  preparation  of  this  report  are  gratefully 
acknowledged. 

Members  of  the  Institute  of  Mathematical  Sciences,  New 
York  University,  who  have  contributed  directly  to  the  project 
are: 

Melvin  P.  Ehrllch,  Assistant  Research  Scientist,  May  I956 

to  Aiogust  1958 
Peter  Kerwin,  Engineer,  February  1956  to  February  1958 
Iren  Livingston,  Technician,  September  195^  to  September 

1958 
John  McNally,  Technician,  September  I958  to  July  1959 
David  M.  Wetstone,  Associate  Research  Scientist,  December 

1956  to  July  1959. 

The  author  wishes  to  acknowledge  with  considerable  apprecia- 
tion the  highly  informative  conversations  he  held  with  Dr. 
Reimar  Liist  (Max  Planck  Institute  for  Physics,  Munich)  which 
contributed  to  parts  of  the  theoretical  discussion. 

The  assistance  of  Mr.  Carl  Bass,  of  the  Institute,  in  the 
preparation  of  many  of  the  diagrams  is  gratefully  acknowledged. 


4  - 


NYO-2877 

TABLE  OF  CONTENTS 

Page 

Al53tract 2 

Personnel  and  Acknowledgments 4 

Table  of  Figures 7 

I.  Introduction  9 

II.  Experimental  Procedures  and  Theory 

A.  General  Description  11 

B.  Creation  of  Plasmoids 12 

C.  Creation  of  Axial  Magnetic  Field 15 

D.  Firing  Synchronization l6 

E.  The  Differential  Solenoid 16 

1.  Normal  Differential  Solenoid l8 

2.  Shorted  Differential  Solenoid  26 

F.  The  Curved  Solenoid 50 

G.  Center-of-Mass  and  Expansion  Velocity  Measurements.  35 

III.  Results  and  Discussion 

A.  General 37 

B.  Plasmoid  Plow  in  the  Normal  Differential  Solenoid  .  37 

1.  Photographic  Results 38 

2.  Photomultiplier  Measurements 46 

3.  Derived  Results 55 

C.  Plasmoid  Flow  in  the  Shorted  Differential  Solenoid.  76 

D.  Plasmoid  Flow  in  the  Curved  Solenoid 77 


-  5  - 


Page 
Appendix  I;   Experimental  Equipment  .  .    87 

Appendix  II;  Coaxial  Plasma  Source.  ...  95 

Appendix  III;  Button  Source  Trigger.  . ,  9^ 

Appendix  IV;  Laboratory  Designed  Circuits  .  97 

Appendix  V;  Photomultlplier  Experiment  .,..,.,.  104 

Appendix  VI;   Determination  of  Flux  Distribution  in  Multiple 
Circular  Coils  with  Axial  Symmetry  ....   107 


-  6  - 


NyO-2877 
TABLE  OF  FIGURES 


Schematic  of  Plasma  Source  with  Differential 
Solenoid  ..... 


o    o    o 


Page 


13 
Ih 

17 
17 
20 


2.  Schematic  of  Plasma  Source  with  Curved  Solenoid. 

3.  Light  Signal  from  Plasma  Source  Mouth.  ..... 

4.  Ringing  of  Plasma  Source  and  Capacitor  Assembly. 

5.  Flux  Line  Pattern  in  Normal  Differential  Solenoid 

6.  Two-Turn  Search  Coll  on  Normal  Differential  Solenoid  .    21 

7.  Fl\«c  Line  Pattern  in  Shorted  Differential  Solenoid  .  .   28 

8.  Two-Turn  Search  Coil  on  Shorted  Differential  Solenoid,    29 

9.  Fl\ix  Line  Pattern  in  Curved  Solenoid  .........    32 

10.  Two-Turn  Search  Coil  on  Curved  Solenoid.  ...  ...  .    33 

11.  Plasmoid  Flow  in  the  Normal  Differential  Solenoid.  I  .   40 

12.  Plasmoid  Plow  in  the  Normal  Differential  Solenoid. II. 

Study  of  Camera  Interpretation  .....  kl 

13.  Plasmoid  Flow  in  the  Normal  Differential  Solenoid. III. 

Effect  of  Field  Variation.  . 43 

14.  Plasmoid  Flow  in  the  Normal  Differential  Solenoid. IV. 

Effect  of  Field  Variation  in  Closeup  ...   44 

15.  Plasmoid  Flow  in  the  Normal  Differential  Solenoid.  V. 

Diverting  of  Plasma  Filament  ...  47 

16.  Plasmoid  Flow  vs.  Flxix  Line  Behavior  at  a  Magnetic 

Mirror  .....,.=.  ..............  .    48 

17.  Photomultiplier  Traces  of  Plasmoid  Transit  without 

PXGXCl*    9    9    «    o    o    «    o    o    •    »    *    e    o    »    •    •    •    •    •    *    •    •    •        T"^ 

18.  Photomultiplier  Traces  of  Plasmoid  Transit  with  Field.   52 

19.  Effectiveness  of  Magnetic  Mirror  for  Variation  of 

Particle  Velocity  Orientation  in  Moving  Plasmoid  .  .   73 


-  7  - 


Page 

20.  Plasmoid  Flow  In  the  Shorted  Differential  Solenoid,  I. 

Single  Jiraper  .  .  .  » 78 

21.  Plasmoid  Flow  in  the  Shorted  Differential  Solenoid.  II, 

Double  Jumper  ..........  .  79 

22.  Plasmoid  Flow  in  the  Curved  Solenoid.  I.   Plasmoid. 

at  4o  Joules 80 

23.  Plasmoid  Flow  in  the  Curved  Solenoid.  II,   Plasmoid  .  , 

24.  Plasmoid  Flow  in  the  Curved  Solenoid.  III.   Plasmoid 

25.  Experimental  Arrangement  with  Differential  Solenoid  .  .  89 

26.  Button  Source  Trigger 90 

27.  Experimental  Arrangement  with  Curved  Solenoid  .....  91 

28.  Schematic  of  System 92 

29.  Detail  of  Coaxial  Plasma  Source  95 

30.  Detail  of  Spark  Gap  Trigger 96 

31.  Pulse  Amplifier  Circuit  .  , ^.  99 

32.  Main  Power  Supply  Circuit  .  , 100 

33.  Ignitron  Unit  Circuit  .....  101 

34.  Control  Unit  Circuit. 102 

35.  Meter  Relay  Circuit IO5 


-  8  - 


NYO-2877 


EXPERIMENTS  ON  SUPERSONIC  PLASMA  FLOW  ALONG  MAGNETIC  FIELDS,  II, 


Introduction 


This  report  presents  the  continuation  of  the  work  first 
described  In  1958  on  the  motion  of  plasmoids  along  diverted 
magnetic  field  lines,  and  is  the  major  and  final  report. 

Recently,  the  experimental  applications  of  plasma  physics 
have  been  largely  the  controlled  thermonuclear  fusion  problem, 
wherein  almost  all  approaches  have  utilized  some  form  of  static 

plasma,  either  self-contained  or  with  a  superimposed  magnetic 

p 
field,  or  both.   The  experiments  reported  here  differ  in  that 

an  individual  burst  of  plasma  (plasmoid)  is  projected  into  an 
evacuated  chamber  and  subjected  to  various  magnetic  field 
geometries.   However,  this  present  work  was  not  directed  toward 
the  attainment  of  fusion  but  was  motivated  primarily  by  the 
desire  to  investigate  a  relatively  new  kind  of  phenomenon  and 
gain  qualitative,  and  some  quantitative,  information  which  might 
serve  to  suggest  and  encourage  research  in  this  direction  of 
plasma  physics. 

Some  fundamental  questions  asked  were,  will  a  uniform  B^ 
field  collimate  plasmoid  flow,  and  over  what  distance?  What 
would  be  the  behavior  of  a  projected  plasmoid  in  regions  of  con- 
verging field  lines  (magnetic  mirror),  diverging  field  lines,  or 
curving  field  lines?  How  do  the  center-of-mass  and  thermal 
velocities  of  the  plasmoid  behave  under  these  conditions?  Some 

of  these  questions  have  already  been  examined  in  the  first 

1 
report,   and  others  are  undertaken  here. 
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"5  4 
Interest  in  plasmolds  is  largely  associated  with  Bostlck  ' 

who  has  described  a  "button"  or  "rail"  source  for  their  genera- 
tion and  propulsion  across  magnetic  fields,  and  has  given 
considerable  thought  to  the  morphology  of  the  plasmold,  itself. 
This  work  was  continued  by  Flnkelsteln,  et  aj^.  ,  who  designed 
and  Initiated  the  presently  reported  program. 

There  are  several  other  investigations  in  which  plasmolds 
were  shot  into  magnetic  fields  or  were  magnetically  accelerated. 
Sinelnlkov,  et^  £l.  ,  employed  magnetic  and  electron  probes  to 
detect  the  radial  polarization  of  a  plasmold  when  entering  an 

axial  magnetic  field,  a  variety  of  plasmold  was  transferred  and 

7 
compressed  in  stages  by  Coensgen  in  a  fusion  device  at  Llvermore 

called  Toy  Top,  and  experiments  were  recently  reported  by 

o 
Marshall,   utilizing  the  technique  of  a  traveling  magnetic  wave 

to  compress  a  plasmold  and  accelerate  it  into  a  vacuum  region. 

To  date,  there  has  been  no  published  attempt  to  guide  a 
plasmold  around  a  curved  field,  as  has  been  tried  here,  though 
there  are  a  number  of  examples  of  sustained,  curved  discharges. 
One  of  the  more  interesting  of  these  is  the  l80°  plasma  beam 
photographed  by  Harrison  and  Dawton-^  at  Harwell,  which  was 
stabilized  by  axial  magnetic  fields. 

It  is  hoped  that  the  broad  survey  of  phenomena ,  investigated 
and  described  so  briefly  herein  and  in  the  previous  report,  will 
serve  to  encourage  further  investigations  in  this  interesting 
subject. 
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II.   Experimental  Procedures  and  Theory 

A.   General  Description;  The  phenomenon  studied  was  the 
behavior  of  a  plasmold  produced  by  a  vacuum  spark  discharge 
which  was  Injected  along  magnetic  field  lines  which  were 
diverted  or  curved,  all  within  an  evacuated  tube.   In  most 
cases  observations  were  made  on  Individual  discharges  (shots), 
although  there  were  indications  that  a  single  shot  may  have  con- 
sisted of  a  number  of  individual  sub-bursts,  due  to  capacitor 
ringing,  which  fused  into  a  single  plasmold  during  passage  along 
the  tube . 

The  plasma  was  created  by  the  high-speed  discharge  of  a 
capacitor  across  a  coaxial  source  (see  B,  below).   The  magnetic 
field  was  created  by  the  discharge  of  a  capacitor  bank  (see  C, 
below)  across  a  suitable  solenoid  around  a  cylindrical  glass 
tube,  the  plasma  being  injected  axially  at  one  end.  The  fields 
reached  maximiom  relatively  slowly  (the  order  of  100  to  I50  (xsec), 
while  the  plasmold  was  created  and  traversed  the  tube  rapidly 
(the  order  of  10  jisec).   Thus,  the  plasmold  saw  an  approximately 
D.  C.  field.   It  was  necessary  to  employ  a  synchronization 
system  to  achieve  this  effect  (see  D,  below).   The  vacuum  equip- 
ment (see  Appendix  I)  maintained  the  pressure  in  the  range  of 
3  X  10~  to  8  X  10"  mm  Hg,  except  during  the  first  few  seconds 
after  each  shot,  but  shots  were  spaced  at  least  one  minute  apart. 

Two  types  of  solenoids  and  tubes  were  employed.   The  first 
was  a  straight  helix  of  constant  diameter,  of  which  the  center 
portion  was  wound  with  smaller  pitch  than  the  ends  (see  B, 
below);  this  helix  was  placed  around  a  straight  cylindrical 
glass  tube.   Most  of  the  observations  reported  herein  employed 
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this  geometry.   A  variation  was  obtained  by  shorting  (and  cutting) 
the  center  section.   The  second  arrangement  employed  a  uniform 
curved  helix  of  constant  diameter,  fitted  around  a  curved  cylin- 
drical glass  tube  (see  F,  below). 

The  bulk  of  the  observations  consisted  of  time-integrated 
photographs  of  individual  plasmoid  shots  traversing  the  tube, 
over  a  wide  range  of  photographic  conditions  and  positions, 
employing  Polaroid  46l  film.   In  addition,  one  geometry  was 
utilized  to  make  an  extensive  study  of  plasmoid  center-of~raass 
and  expansion  velocities,  employing  the  photomultiplier  technique 
(see  G,  below). 

As  will  be  noted,  the  details  of  the  experimental  equipment 
will  be  found  in  various  parts  of  the  Appendices. 

B.  Creation  of  Plasmoids;  The  plasmolds  were  created  by  a 
vacuum  spark  discharge  across  a  coaxial  source  used  in  conjunc- 
tion with  a  low-inductance  capacitor.  Mechanical  details  of  the 
system  will  be  found  in  Appendix  II.  The  capacitor  was  rated  at 
0.55  M.P  and  0.040  ^h,  and  was  fired  at  8,  12,  and  I6  kV. 
However,  unless  otherwise  stated,  firing  was  at  12  kV,  producing 
about  40  Joules. 

The  plasma  source  was  positioned  coaxially  at  one  end  of 
the  tube  and  solenoid,  the  positioning  being  somewhat  different 
for  each  of  the  two  solenoid  geometries  employed.   These  two 
arrangements  are  shown  schematically  In  Figures  1  and  2.   After 
some  Initial  trials  it  was  deemed  necessary  to  surround  the 
source  with  a  short  glass  tube  for  purposes  of  initial  collimation 
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(a  technique  subsequently  reported  also  by  Sinelnlkov,  £t  a^i'  )» 
this  in  no  way  vitiated  the  genuineness  of  the  colllmating 
effect  provided  by  the  magnetic  field,  as  is  amply  demonstrated 
by  photographs  in  which  the  solenoid  was  not  energized. 

It  was  of  some  interest  to  investigate  the  time  pattern  of 
discharge  at  the  source  and  to  compare  it  to  the  ringing  of  the 
assembly.   Examples  of  these  phenomena  are  shown  in  Figures  5 
and  4.   In  order  to  eliminate  the  effects  of  noise,  the  source 
picture  was  taken  as  the  output  of  a  photomultiplier  at  the  end 
of  a  twenty  foot  acrylic  light  pipe  which  was  pointed  directly 
at  the  mouth  of  the  source.   The  frequency  measured  was  1.29  mc. 
The  ringing  frequency  of  the  assembly  was  measured  directly 
across  the  capacitor j  it  was  620  kc.   (The  imperfect  wave  shape 
is  probably  due  to  poor  matching  of  the  voltage  attenuation 
chain. )  This  is  in  satisfactory  agreement  with  the  assumption 
that  the  source  fired  a  sub-burst  for  a  pair  of  successive  half- 
cycles  or  full  cycles  of  the  assembly  (each  light  peak  being 
not  necessarily  a  plasma  sub-burst).   A  choice  between  these 
possibilities  is  afforded  by  the  data  of  Section  IIIB2. 

The  brass  electrodes  of  the  plasma  source  were  assximed  to 
produce  a  copper-zinc  plasmoid,  with  all  the  ions  taken  at 
64  g/g-atom  for  purposes  of  computation. 

C.   Creation  of  the  Axial  Magnetic  Field;  The  field  was 
created  by  one  of  several  solenoid  geometries  (see  E  and  F, 
below),  within  the  experimental  tube.   The  pulsed  current  re- 
quired was  obtained  by  the  discharge  of  a  capacitor  bank  across 
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the  solenoid  using  Ignltron  switching  (for  mechanical  and 
electrical  details,  see  Appendices  I  and  IV).   No  special 
attempt  was  made  to  reduce  external  inductance  as  high  speed 
was  not  a  requisite. 

The  bank  consisted  of  eight  105  M-F  oil  and  paper  capaci- 
tors of  approximately  0.01  ^.h  each,  in  parallel  and  normally- 
charged  to  3  kV^  total  energy  was  therefore  3780  Joules.   They 
were  connected  by  a  pair  of  aluminum  sheets  separated  by  Mylar. 
(Leakage  resistance  was  600  megohms. )  The  external  inductance 
(total  circuit  less  bank)  was  of  the  order  of  8  to  10  iih 
(depending  on  solenoid  geometry)  with  about  two-thirds  contri- 
buted by  the  solenoid. 

D.  Firing  Synchronization;   A  two-turn  search  coil  was 
placed  around  each  solenoid.   Initial  rise  in  this  coil  was 
extremely  fast  and  was  taken  as  "time  zero"  for  purposes  of 
synchronization!  return  to  zero  amplitude  represented  attain- 
ment of  maximum  field  in  the  solenoid.   The  time-zero  signal 
was  used  to  initiate  a  time-delay  mechanism  (see  Appendices  I 
and  rv)  for  firing  the  plasma  source,  so  arranged  that  the 
plasmold  reached  approximately  the  center  of  the  solenoid  at 
maximum  field.   Jitter  in  this  circuit  was  permissible  within 
wide  latitude  because  the  field  at  maximum  changed  slowly  (first 
derivative  very  small)  during  the  plasmold  transit  time. 

E.  The  Differential  Solenoid:  This  coil  was  constructed 


of  0.12  inch  diameter,  "ha If -hard"  copper  wire,  wrapped  tightly 
around  a  50-lnch-long  cylindrical  Pyrex  pipe,  2.0  inch  I.D., 
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Figure  3 

Light  Signal  from  Plasma 
Source  Mouth 
(  1  txsec/cm  ) 


Figure  ij. 

Ringing   of   Plasma  Source 
and   Capacitor  Assembly 
(    1  tisec/cm   ) 
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2.5  inch  O.D. ,  producing  a  current  sheath  of  1.31  inch  radius. 
The  first  seven  turns  were  spaced  at  lolO  inches  per  turn,  the 
next  nine  turns  at  O.hj)   inch,  and  the  remaining  eleven  turns 
again  at  1.10  inches.   The  coll  was  utilized  in  two  different 
configurations,  which  will  be  discussed  separately. 

•  1.   Normal  Differential  Solenoid;  The  bulk  of  the 
observations  were  obtained  with  the  solenoid  employed  as  des- 
cribed above,  namely  with  a  central  region  in  which  the  field 
intensity  was  more  than  double  that  of  the  two  end  regions. 
The  flux  line  pattern  is  indicated  roughly  by  Figure  5-   This 
provided  the  plasmoid  with  an  entrance  into  and  an  exit  from  a 
magnetic  mirror  during  its  transit.   The  search-coll  time  for 
this  arrangement  is  shown  in  Figure  6  (amplitude  Is  in  arbitrary 
units  and  is  not  the  same  for  both  pictures).   From  6a.  one 
obtains  132  fisec  as  the  time  to  reach  maximum  field  from  time 
zero.   The  decay  constant,  a,  for  this  damped  sinusoidal  wave 
may  be  obtained  from  6b, ,  utilizing  the  data  of  the  first  two 
extrema  °. 

tn  =  0  ^2  "^  ^^^  ^^®^ 

v,    =  77  units  Vp  =  60  units 

«  =  2M30   ^"^  S  =  9.^4  X  10^  sec"^ 

This  Is  also  the  decay  constant  for  the  solenoid^  the  measured 
frequency  is  1.90  kc. 

It  is  of  importance  to  determine  the  field  strength  (i.e., 
magnitude  of  magnetic  induction)  within  the  solenoid.   It  was 
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not  convenient  to  utilize  the  direct  approach  of  placing  a  small 
search  coil  within  the  experiment  tube.   The  simplest  alternati-ve, 
therefore,  would  be  to  integrate  the  outside  search-coil  trace 
(Jvdt),  yielding  total  contained  fl\ix,  and  to  devide  by  the 
current  sheath  cross-section.   This  suffers  from  two  disadvan- 
tages.  The  first  is  that  ascribing  to  the  solenoid  the  flux  seen 
by  the  search  coil  in  effect  assiomes  unit  coupling  between  the 
two  coils.   Some  careful  measurements,  corroborated  by  theoreti- 
cal calculations  of  mutual  inductance,  indicated  values  of  the 
coupling  corfficient  below  0.3,  making  this  procedure  highly 
uncertain.   Second,  the  distribution  of  flux  is  not  uniform  with- 
in the  solenoid  J  although  these  helices  are  not  bad  approximations 
to  true  current  sheaths,  one  is  not  certain,  a  priori,  how  close 
the  assumption  would  be.   Thus,  even  if  total  flux  were  obtain- 
able, division  by  the  cross-section  would  produce  only  an 
average  field  intensity  of  uncertain  validity. 

An  alternate  method  would  be  to  calculate  the  self-inductance 
which,  combined  with  a  knowledge  of  current,  would  yield  values 
for  the  total  flux.   Calculation  of  field  intensity  would  still 
suffer  from  the  second  objection  noted  above.   Calculation  of 
inductance  for  a  coil  geometry  of  this  sort  is  difficult.   It 
can  be  done  readily  for  the  three  parts  of  the  solenoid  on  an 
individual  basis,  but  the  mutual  end  effects  are  definitely  not 
second  order. 

In  any  case,  a  knowledge  of  the  maximtim  current  would  be 
required.   For  this  solenoid  a  direct  measurement  was  not  made, 
but  the  circuit  parameters  obtained  by  the  search  coil,  plus 
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Figure   6 
Two-Turn  Search  Coll   on  Normal  Differential  Solenoid 
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bank  capacitance  and  initial  voltage,  are  sufficient  to  yield  a 
reliable  value  by  calculation.  The  general  circuit  equation  in 
this  case  is 

Boundary  conditions  are  that  1=0  for  t  =  0  or  infinity.   For 
an  under-damped  (ringing)  circuit,  the  solution  is 


/l    r2  ^  1/2 
where  a  =  7^  and  a,   =  2irf   -  ,%  -  ^^ '    •   (Since  a  and  f  have 

been  measured,  and  C  is  known,  one  may  calculate  L  and  Rs 


(27rf )2  -^-  o^ 


L  - 


C[(2Trf)^  +  a^] 


8.4  X  10"^[(2Tr  1.9  X  10^)^  +  (9.^^  X  10^)^] 


=  8.3  jih 


and 


R  =  2La  =  2(8.3  X  10"^)(9=^4  X  10^) 


=  1.57  X  10"^  ohm. ) 


It  may  be  reasonably  assumed  that  the  bank  is  fully  dis- 
charged when  the  current  has  reached  maxlmiom  at  t^.   Thus 
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The  value  of  the  function  within  the  brackets,  for  t,  =  132  M.sec, 
is  8.92  X  10"^.   if  the  initial  voltage  on  the  bank  is  3  kV,  then 

1    ,    ^o      ^  (8.4  X  10~^)(3  X  10^) 
'"^'^  '  8.92  X  10"^       8.92  X  10-5 

=  28,300  amp. 

This  is  the  undamped  maximvmi  current;   the  actvial  current,   i,  ,   at 
t,    is  given  by 

-at, 

i,    =  i  e        ^ 

1         max 

=  24,900  amp. 

As  an  alternative  to  a  direct  calculation  of  self-inductance 
(which,  as  noted  above,  would  be  not  only  difficult  but  of 
limited  utility)  the  following  procedure  wa^  adopted.   The 
solenoid  was  approximated  by  a  model  consisting  of  coaxial  cir- 
cular current  filaments  of  equal  radii,  spaced  at  intervals  which 
corresponded  to  the  actual  turns  of  the  solenoid,  twenty-seven 
filaments  in  all.   Since  this  resulted  in  cylindrical  syTmnetry 
the  problem  could  be  reduced  to  two  dimensions  by  projection  on 
the  r-z  plane,  with  further  reduction  to  the  r  >  0  half.   Equations 
were  derived  from  the  Biot-Savart  law  for  the  longitudinal  and 
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radial  components  of  the  magnetic  Induction,  B  and  B  ,  for  each 
filament  at  any  point  on  the  r-z  plane  (and  by  symmetry,  for  any 
point  in  space),  and  then  summed  over  all  filaments: 

IX  i    2r  (r  -  r) 
a 

M.  i   z-z    2rr 
a 

where  (r,z)  is  the  point  in  question,  (r  ,z  )  the  intersection 
of  the  filament  in  the  plane,  and  K  and  E  are  the  complete 
elliptic  integrals  of  the  first  and  second  kind,  respectively. 
Also  one  has 

A^  =  (r  +  r^)  +  (z  -  Zq)   , 

a     =  (r  -  r^)   +  (z  -  z^)   , 

2   ^^^o 
and      k  =  — 5— 

X 
was  the  modulus  of  the  integrals.   The  problem  was  Fortran  coded 
for  the  IBM  704  computer  from  which  were  obtained  the  values  of 
B  at  various  significant  positions  on  the  half -plane  and  thus, 
by  symmetry,  within  the  coil  (see  Appendix  VI),   The  pattern  of 
flux  distribution  was  obtained  also,  from  which  Figure  5  was 
drawn.   The  problem  was  run  for  unit  current. 

Since  the  model  was  based  on  filaments  carrying  unit  current 
normal  to  the  solenoid  axis,  the  values  of  B  obtained  should  be 
multiplied  by  the  radial  component  of  the  projection  of  the 
helical  current  on  the  r-z  plane  (i.e.,  the  aximuthal  current). 
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This  will  differ  slightly  for  region  b.  as  compared  to  a.  and 
c.  (Figure  5).   Since  the  distance  for  a  half  turn  (regions  a. 
and  c.)  was  0.55  inch  while  the  diameter  was  2.62  inches,  the 
szimuthal  component  of  the  current  is  given  by 

i^  =  i  cos  tan"-"-  S^ 
az  2.62 

=  0.98  i  . 

For  a  helical  current  of  24,900  amp,  the  azimuthal  component 
would  be  24,400  amp.  In  region  b.  the  correction  amounts  to 
less  than  half  a  percent,  and  i   =  i,  =  24,900  amp. 

This  result  may  now  be  combined  with  that  of  the  machine 
computation  to  yield  field  strengths.  On  the  axis  the  magni- 
tude of  magnetic  induction  was  approximately; 

central  portions  of  regions  a.  and  c. ;   1.1  weber/m 
central  portion  of  region  b.         1   2.5  weber/m  . 
Attempts  to  simplify  the  picture  of  the  field  distribution  off 
the  axis  suffer  difficulties  whether  one  regards  the  model  or 
the  actual  helical  coil.   In  the  former  case  moving  radially 
from  the  axis  results  in  either  a  decrease  or  an  Increase  in 
field  intensity  depending  upon  whether  z  is  located  between  two 
circular  filaments  or  within  one.   In  the  latter  case  (considering 
the  helix  itself),  the  variation  in  intensity  is  dependent  upon 
the  angular  direction  of  r  for  a  given  z.   Thus  the  model 
employed  here  becomes  increasingly  invalid  at  substantial  dis- 
tances from  the  axis.   However,  values  of  B,  obtained  from  the 
machine  data,  undergo  variations  of  no  more  than  5  percent  at 
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r  =  0.5  inch  for  the  central  positions  of  the  a.  and  c.  regions 
and  2  percent  for  the  central  portion  of  the  b.  region.   Most  of 
the  plasmold  flow  lay  well  within  this  radius. 

2.   Shorted  Differential  Solenoid;  It  was  of  particu- 
lar interest  to  obtain  some  indication  of  the  behavior  of  the 
plasmold  as  it  emerged  from  and  reentered  a  magnetic  mirror 
(flow  through  a  magnetic  bottle).  This  was  accomplished  by 
shorting  (and  cutting)  region  b.  of  Figure  5.   The  flux  line 
pattern  is  indicated  roughly  by  Figure  7-  Two  shorting  Jumpers 
are  Indicated i  some  photographs  were  taken  with  one  connected, 
some  with  both,  as  will  be  noted.   The  search  coil  trace  for 
this  arrangement  is  shown  in  Figure  8  (again,  amplitude  is  in 
arbitrary  units  and  is  not  the  same  for  8a,  and  8b.).   From  8a. 
one  obtains  108  iisec  as  the  time  to  reach  maximum  field.  The 
decay  constant  may  be  calculated  from  8b. ,  utilizing  the  data 
of  the  second  and  third  extrema  (the  first  being  off-scale); 

tg  =  200  iisec  t,  =  4l6  usee 

Vg  =  87  units  V,  =  57  units 

"  "  416-^00  -«n  ^  =  1.21  X  10-^  sec"  . 

The  corresponding  frequency  is  2.31  kc. 

From  these  results,  the  determination  of  L  and  R  follow 
(as  described  in  Section  El,  above);  L  =  5.6  [ih; 
R  =  1.36  X  10   ohm.   Calculation  of  current,  as  in  El,  yields 
^mo-ir  -   3^,200  amp,  from  which  the  actual  current  at  maximum 
field  is  given  by  1^  =  30,000  amp  at  t^  =  108  p-sec. 

This  problem  also  was  submitted  to  computer  analysis  in  a 
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fashion  similar  to  that  of  the  normal  Differential  Solenoid. 
Eight  circular  filaments  represented  region  a.  and  twelve, 
region  c,  with  an  appropriate  space  between  (the  now  inopera- 
tive turns  of  region  b.)-  The  pattern  of  flux  distribution  is 
that  shown  roughly  in  Figure  7.  The  values  of  B  obtained  were 
again  for  unit  current  and  were  multiplied  by  the  azimuthal 
component  of  1,,  above,  namely  1,  reduced  by  two  percent? 
^az  ^   29,400  arap.   On  the  axis  the  magnitudes  of  magnetic 
induction  were  approximately? 

central  portions  of  regions  a,  and  c.t   1=5  weber/m 
central  portion  of  region  b.         s  0„4l  weber/m  , 
Variations  of  B  off  the  axis  for  the  central  portions  of  regions 
a.  and  c.  were  less  than  3   percent  for  r  =  0.5  inch  or  less. 

The  situation  in  region  bo  is  somewhat  more  complicated,  in 
part  because  the  plasmoids  expanded  to  the  wall  of  the  tube 
(r  =  1.0  inch)  and  in  part  because  a  curving  transverse  field 
of  some  magnitude  was  created  by  the  Jxorapers  (see  Figure  7).   On 
the  axis  at  the  center  of  the  b,  region  the  flxix  density  was  0.4l 
weber/m  ,  as  noted  above,  decreasing,  in  the  radial  direction,  by 
6.8  percent  at  r  =  0.5  inch  and  by  25  percent  at  the  wall  (r  =  1.0 
inch).   Furthermore,  the  jumpers  created  azimuthal  fields,  normal 
to  the  plasmoid  flow,  of  approximately  0.25  weber/m  at  r  =  1.0 
for  one  jumper  and  half  that  value  when  two  jumpers  were  employed! 
these  values  decreased,  of  course,  proceeding  inward  toward  the 
axis.   This  is  to  be  compared  to  0.51  weber/m  at  the  inner  wall 
(r  =  1,0),  resulting  from  the  solenoid  current. 

One  additional  set  of  calculations  is  important.   For  the 
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Figure   8 
Two-Turn  Search  Coil   on  Shorted  Differential  Solenoid 
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Curved  Solenoid  (to  be  discussed  below)  there  was  no  syminetrlc 
model  available  susceptible  of  straightforward  analysis.  Thus, 
it  was  necessary  to  obtain  a  rough  measure  of  field  strength, 

i.e.,  flux  density,  through  calculation  of  coil  inductance.   It 
is  of  interest  here  to  ascertain  the  degree  of  reliability  for 
such  a  procedure  in  a  case  where  independent  data  on  flux 
density  can  be  obtained i  specifically,  the  last  ten  tiirns  of 
region  c.  in  the  Differential  Solenoid  were  treated  as  an 
Independent  coil  and  the  inductance  was  measured  directly  by  a 
determination  of  reactances  L  -  I.56  ^JLh.  Therefore,  total 
flux  is  given  by 

X   ^az^  „  (2.9^  X  10^)(1.56  X  10"^) 

-  4.59  X  10"-^  weber. 

Cross-section  for  r  =  I.31  inch  is  3-50  x  10  ^  m  .   Therefore, 
B   =  4.59  X  10"V3»5  X  10"^  =  1.31  weber/m^o   This  is  to  be 
compared  to  a  machine  calculation  of  axial  field  strength  for  an 
eleven  circular  filament  model,  which  yielded  B  -   I.50  weber/m  , 
The  correspondence  is  quite  fortuitous  1  there  are  a  number  of 
errors  inherent  to  this  procedure  but  fortunately  many  of  them 
cancel,  suggesting  that  it  is  feasible, 

F.   The  Curved  Solenoids  The  other  major  geometry  investi- 
gated was  a  solenoid  which  created  a  curved  field.   The  Pyrex 
experiment  tube  was  about  45  inches  long,  with  2  inch  I.D.  and 
2.3  inch  O.D.   Thirty-one  turns  of  0.12  inch  diameter,  "half- 
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hard"  copper  were  wrapped  tightly  around  at  a  mean  arc  spacing 
of  1.05  Inches  per  turn,  producing  a  current  sheath  of  1.20  inch 
radius.  This  solenoid  approximated  a  6o°  arci  radius  of  curvature 
for  this  geometry  (as  well  as  the  length  of  the  chord)  was  31 
inches.  The  flvix  pattern  is  Indicated  roughly  in  Figure  9.   The 
field  is  not  distributed  with  axial  symmetry,  of  course,  being 
somewhat  more  concentrated  toward  the  center  of  curvature;  how- 
ever, the  effect  is  not  as  exaggerated  as  in  the  figure.   For  a 
rough  estimation  of  field  distribution  in  geometries  of  this 
sort  one  may  assume  that  the  product  rB  is  a  constant,  where  r 
is  the  distance  from  the  center  of  curvature  to  some  point  within 
the  coil  on  the  curvature  plane.   From  this  one  may  readily  deduce 
that  the  logarithmic  gradient  of  B  (i.e.,  |V!b|/b)  is  roughly  equal 
to  the  reciprocal  of  r  at  that  point.   Thus,  variation  in  flux 
density  due  to  ciorvature  should  amount  to  approximately  6  percent 
within  the  tube. 

The  search-coil  trace  for  this  arrangement  is  shown  in 
Figure  10  (again,  amplitude  is  in  arbitrary  units  and  Is  not  the 
same  for  10a.  and  10b.).   Figure  10a.  shows  the  time  to  reach 
maximum  field  as  120  M-sec.   The  decay  constant  is  again  calculated, 
from  Figure  10b,  using  the  first  and  second  extrema s 

t,  =  0  t^  =  240  M-sec 

v^  -  80  units  Vp  =  57  units 

a  -  ^^   -«n  |y  =  1.41  X  10^  sec"^. 

The  corresponding  frequency  is  2.09  kc.   Values  of  L,  R,  and 
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Figure   10 
Two-Turn  Search  Coil   on  Curved   Solenoid 
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max 


are  again  calculated  as  In  Section  El,  above  i  L  =  6.8  \ih; 


_2 
R  =  1.9  X  10   ohmj  l^g^  =  50,800  amp.   The  actual  current  at 

maximum  field  was  i^  =  26,000  amp  at  t^  =  120  ^JLsec.   Azimuthal 
current  Is  again  about  two  percent  less,  or  i   =  25,400  amp. 
It  Is  of  interest  to  note  that  in  this  instance  an  independent 
measure  of  current  was  obtained  to  determine  the  validity  of 
these  mathematical  procedures.   A  strip  of  manganin  metal  was 
placed  in  series  with  this  solenoid  and  the  potential  drop 
recorded  by  an  oscilloscope  during  bank  discharge.   At  3  kV,  th« 
maximum  of  the  trace  was  29.0  V.  The  resistance  of  the  manganin 
strip  was  measured  independently,  yielding  1.11  x  10"^  ohm. 
The  current  at  maximum  discharge  was  therefore  26,100  amp,  which 
corresponds  quite  closely  to  i^  above. 

As  has  been  discussed  in  Section  E2,  above,  there  is  no 
model  for  this  solenoid  which  possesses  sufficient  symmetry  to 
permit  straightforward  determination  of  flux  distribution  even 
by  machine  computation,  at  least  with  the  degree  of  relative 
simplicity  permitted  by  the  Differential  Solenoid.   It  was  there- 
fore necessary  to  proceed  by  the  approximation  method  already 
discussed,  and  shown  to  be  reasonably  valid,  employing  measured 
inductance.   From  an  inductive  reactance  measurement  on  this 
coil,  L  =  4.87  ijLh  for  J>1   turns.   Therefore,  total  flux  is  given  by 

¥  =  ^^^^  =  (2-54  X  10^H4.87  X  10"^) 

=  3.99  X  10"-^  weber. 

Cross-section  for  r  =  1.20  inch  is  2.92  x  10~^  m^.   Therefore, 

^av  "  ^°95  ^  10"V2.92  X  10"^  =  1.37  or  approximately  1.4  weber/m^. 
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G.   Center-of-Mass  and  Expansion  Velocity  Mea svirement s  s 
Two  parameters  of  the  plasmold  which  could  be  obtained  In  a  fairly 
straightforward  fashion  were  velocity  of  the  plasmold  center-of 
mass  and  the  rate  of  expansion  of  the  plasmold  during  flight. 
Both  weare  obtained  as  the  result  of  the  same  series  of  measure- 
ments, namely  tlme-of -flight  studies  with  photomultlpllers,  using 
the  normal  Differential  Solenoid, 

Preliminary  observations  suggested  considerable  variability 
in  timing  with  respect  to  various  attempted  methods  of  osclllo- 
iscope  triggering  ("Jitter"),  To  eliminate  the  trigger  as  a  time 
marker,  three  photomultlpllers  were  employed,  one  focused  at  about 
the  third  turn  of  region  a,,  8  cm  from  the  mouth  of  the  source, 
the  second  at  the  boundary  of  regions  a.  and  b.,  and  the  third 
at  the  boundary  of  regions  b,  and  c.   The  resulting  distances  of 
plasmold  traverse  were  11. 5  cm  and  10.5  cm  for  regions  a,  and  b., 
respectively.  Thus,  two  different  sets  of  results  were  obtained 
sj).multaneously  for  each  plasmold  burst.  The  photomultlpller  out- 
puts were  run  into  three  oscilloscopes  of  the  same  model  and  make, 
the  sweep  generators  of  which  had  been  mutually  calibrated,  and 
suitable  triggering  was  employed  to  minimize  trace-time  error 
(see  Appendix  V  for  further  details). 

Calculation  of  center-of-mass  velocity  depended  upon  the 
assumption  that  the  center  of  mase  would  coincide  at  all  times 
with  location  of  maximum  recombination  light  intensity.   Expansion 
velocity,  on  which  maximvim  temperature  was  based,  was  deduced  from 
the  rate  of  longitudinal  expansion  of  the  photomultlpller  trace 
between  each  of  the  two  pairs  of  locations.   Usually  this  was 
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calculated  from  the  increase  in  half-width  of  a  light  peak, 
measured  at  half-amplitude  (to  be  discussed  further  with  the 
data  in  Section  III).  While  this  method  of  ion-temperature 
determination  is  not  trustworthy,  compared  to  presently  avail- 
able methods,  it  gains  considerably  in  simplicity,  particularly 
where  only  an  upper  bound  is  required. 
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III.  Results  and  Discussion 

A.  General;  A  large  portion  of  the  reported  data  consists 
of  time-Integrated  photographs  of  plasmold  behavior  when  they  are 
injected  in  the  normal  Differential  Solenoid.  The  pattern  of 
the  flow  is  compared  to  that  of  the  field  configuration.  Utiliz- 
ing this  same  arrangement,  tlme-of -flight  data  were  obtained, 
from  which  variations  in  velocity  and  temperature  are  deduced.   A 
method  of  shock-heating  the  original  spark  channel  is  suggested. 
On  the  basis  of  these  results,  and  some  additional  measurements, 
a  value  for  density  is  obtained.   Assiimptions  of  electron  tempera- 
ture are  siiggested,  from  which  a  degree  of  ionization  is  calculated. 
From  this,  Larmor  radius  and  drift  velocity  follow.   A  total 
reflection  coefficient  for  the  plasmold  is  calculated,  viewing  the 
constriction  of  field  in  the  b.  region  as  a  magnetic  mirror.   In 
addition,  an  estimation  is  made  of  the  magnitude  of  the  plasmold 
self-field,  as  compared  to  the  superimposed  B  field. 

Time-integrated  photographs  are  then  shown  of  the  shorted 
Differential  Solenoid  under  certain  varied  conditions.   Finally, 
photographs  are  shown  of  plasmold  behavior  in  the  Curved  Solenoid 
with  variations  of  both  plasmold  and  field  energy,  and  the  effects 
are  calculated  of  an  internal  polarization  field  due  to  charge 
separation. 

B.   Plasmold  Flow  in  the  Normal  Differential  Solenoid;  This 
section  is  divided  into  photographic  results,  photomultlplier 
measurements,  and  derived  results.   The  axial  values  of  magnetic 
induction  are  given  in  each  case,  for  regions  a,  and  c.  and  for 
region  b.  These  values  were  varied  somewhat  for  the  photographic 
results  (by  changing  initial  voltage)  to  accentuate  certain 
phenomena.  The  plasmold  source  was  fired  at  12  kV,  producing  an 
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energy  per  burst  of  40  Joules,  Synchronization  was  arranged  to 
place  the  burst  approximately  10  cm  within  region  a.  at  maximum 
field. 

1.   Photographic  Results;  Pictures  were  taken  of 
Individual  bursts,  with  shutter  open  and  aperture  adjusted  to 
minimize  peripheral  light  (slower  fragments  of  the  plasmoid  that 
had  expanded  to  the  walls)  and  accentuate  the  flow  path  of  the 
central  portion  of  the  burst.  The  photographs  shown  here  are 
representative  of  an  extensive  amount  of  reproducible  data. 

An  example  of  the  flow  pattern  of  the  plasmoid  is  given  in 
Figure  11,  in  which  are  included  a  camera  field-of-view  picture 
of  the  solenoid  and  a  picture  of  plasmoid  discharge  without 
magnetic  field. 

In  this  type  of  photography  the  appearance  6f  the  plasmoid 
flow  in  the  picture  (for  any  given  set  of  experimental  conditions) 
is  strongly  influenced  by  camera  aperture,   A  series  of  pictures 
was  taken  with  gradually  decreasing  aperture  while  holding  all 
other  conditions  constant.  These  are  shown  in  Figure  12.  Not 
only  the  brightness  but  the  curvature  of  the  flow  outline  is 
altered  on  diminishing  aperture  (Figures  12g.  and  12h, )  The 
central,  denser  portions  of  the  plasmoid  pass  into  the  mirror 
between  regions  a.  and  b.  with  little  constriction.  The  fact 
that  the  photographs  are  time  integrated  raises  the  possibility 
(further  suggested  by  time-of=f light  studies,  below)  that  the 
leading  edge  or  first  sub-burst  of  the  plasmoid  is  considerably 
less  affected  by  the  mirror  than  the  remainder,  and  it  may  be  the 
passage  of  this  portion  which  is  being  recorded  in  the  smaller 
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aperture  pictures » 

A  datum  which  will  be  of  later  interest  may  also  be 
obtained  from  this  series  of  figures,  namely  an  estimate  of 
the  plasma  drift  across  a  uniform  field.   Figure  12f.  is 
most  convenient  for  this  purpose  since  it  appears  to  show 
largely  the  flow  of  the  leading  edge  or  high  density  sub- 
burst  (as  described  above).  There  are  six  turns  visible 
in  region  Co  of  this  figure.   If  one  examines  the  flow 
between  the  second  and  sixth  turns  (where  the  field  is 
reasonably  uniform)  one  finds  an  increase  in  plasmoid  diameter 
of  approximately  two  percent,  or  a  transverse  drift  velocity 
of  about  one  percent  of  the  axial  velocity. 

The  next  point  of  interest  concerns  variations  in  flow 
pattern  when  the  field  is  changed.   Figure  1J>   shows  an 
example  of  the  effect  of  doubling  the  field.   At  the  lower 
field  strength  there  is  a  constriction  in  region  b,|  while 
doubling  the  field  produces  a  smaller  beam  in  region  a,,  it 
does  not  appear  to  constrict  further  the  dense  portion  of  the 
plasmoid  flow  in  region  b.   It  does,  however,  reduce  the  hazy 
light  (presumably  the  cooler,  trailing  portions  of  the  plas- 
moid).  Apparently  there  is  a  limit  to  the  degree  of 
constriction  obtainable  under  these  particular  conditions. 
This  phenomenon  was  examined  with  close-up  photography,  as 
is  shown  in  Figure  14.  Here  (in  l4d, )  it  is  seen  clearly  that 
there  is  a  dense  portion  of  the  plasmoid  which  appears  to 
experience  no  further  constriction  when  entering  region  b,,  if 
the  field  is  already  beyond  some  (unspecified)  magnitude. 
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Figure  11 
Plasmold  Flow  in  the  Normal  Differential  Solenoid. I, 
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Plasmoid  Flow  in  the  Normal  Differential  Solenoid, 
II,   Study  of  Camera  Interpretation 
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Figure  12  continued 
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Figure  13 

Plasmold  Flow  in  the  Normal  Differential  Solenoid. 
III.   Effect  of  Field  Variation 


-  i|3  - 


Solenoid 


li^.a, 


Tmi 


««•  IMI 


li^b. 
Plasma:  14.0  joules 
No  field 


ll;c. 
Plasma:  UO  joules 
Regions  a,  and  c: 

0,55  weber/m2        _ 
Region  b.:  1,25  webers/m 
Aperture:  f  16  -  22 


««•  Ilk 


ll4.d. 
Plasma:  i|0  joules 
Regions  a,  and  c: 

1.1  webers/m2 
Region  b,:  2.5  webers/m 
Aperture:  f  22 


Figure  llj. 

Plasmoid  Flow  in  the  Normal  Differential  Solenoid. 
IV.   Effect  of  Field  Variation  in  Close-up 
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Another  feature  of  Figure  l4d.  Is  the  stria t ion  of  the   | 
plasmold  flow  .    A  filament  of  plasma  across  the  top  of  the 
tube  Is  diverted  Inward  upon  entering  the  mirror.  This  appears 
more  clearly  In  Figure  I5,   It  is  not  easy  to  decide  whether 
this  filament  Is  part  of  the  Initial,  high  density  sub-burst 
(which  Is  seen  more  clearly  In  Figure  Ikd, )   or  Is  part  of  the 
trailing  material  which  makes  up  the  general  haze  surrounding 
the  main  beam,  although  the  latter  is  more  likely. 

The  data  from  machine  computation  described  in  Section  IIEl, 
above,  provide  an  opportunity  to  examine  some  details  of  the 
plasmold  flow  in  comparison  to  the  known  behavior  of  flux  lines 
(assuming  accuracy  of  the  model  employed).   Figure  l4c,  provl'des 
a  reasonably  clear  picture  of  plasmold  constriction  upon  entering 
the  mirror  region^  Figure  15  is  particularly  interesting  for  the 
presence  of  the  single  filament  of  plasma.  These  two  flow 
patterns  have  been  transposed  to  Figure  I6,  together  with  points 
(marked  x)  tracing  three  flux  lines  which  lie  in  the  immediate 
vicinity  of  the  pattern  outlines,   A  fl\ix  line  (actually,  the 
trace  of  a  stream  function)  is  to  be  interpreted  as  the  generator 
of  a  surface  of  axial  revolution  upon  which  there  is  no  normal 
component  of  B  (l,eo,  no  flxix  crosses  this  surface).   Radial 
coordinates  are  relative  to  the  center  of  the  solenoid  and  axial 
coordinates  are  relative  to  the  position  of  the  first  turn.  The 
locations  of  turns  (per  the  model)  are  shown.  Magnitudes  of  B 
are  radially  constant  to  better  than  one  percent  within  the 
volxjme  occupied  by  the  plasmold.   Actual  values  of  flux  and  mag- 
netic induction  are  not  given  on  the  diagram  as  it  is  composed 
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of  data  taken  at  two  different  values  of  solenoid  current. 

The  correspondence  between  the  outer  flux  line  and  the 
plasma  filament  is  good,  that  between  the  center  and  inner  lines 
and  the  outlines  of  the  central  flow  of  the  plasmoid  is  not, 
particularly  to  the  left  where  there  is  a  greater  uncertainity 
in  the  location  of  the  plasmoid  flow  boundary.   Some  of  the 
discrepancy  may  be  laid  to  difficulties  of  transposition  from 
photograph  to  figure  and  some  to  the  small  scale  errors  due  to 
employing  the  circular  current  filament  model.   But  primarily 
the  problem  involves  a  decision  as  to  the  location  of  the  bound- 
ary in  region  a,  by  which  we  might  define  the  outline  of  the 
central,  dense  portion  of  the  flow.   One  may  estimate  its  radius 
to  be  about  0,2  inch  or  more  in  region  a,  and  about  0,12  inch  in 
region  b, 

2,   Photomultiplier  Measurements;  The  techniques 
employed  were  described  in  Section  IIG,   All  data  reported  here 
were  of  a  single  40  Joule  plasmoid  fired  into  the  tube  without 
field  or  with  a  field  of  1.1  webers/m  on  the  axis  of  the  a.  and 
c.  regions  and  2,5  webers/m  on  the  axis  at  the  center  of  the  b. 
region.   A  representative  series  of  three  oscilloscope  traces  is 
shown  in  Figure  17  for  a  plasmoid  transit  without  field  (signals 
are  negative-going  and  of  arbitrary  amplitude).   As  previously 
stated,  these  represent  points  located  toward  the  beginning  of 
region  a.,  at  the  a.-b,  transition  and  at  the  b.-c,  transition, 
with  separations  of  11.5  cm  and  10.5  cm,  respectively. 

Figure  17a.  suggests  the  presence  of  a  double  peak  in  the 
plasmoid  but  as  it  progresses  down  the  tube  this  distinction 
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Figure  17 
Photomultiplier  Traces  of  Plasmoid  Transit  without  Field 
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vanishes  and  the  plasmold  becomes  continuous.   This  was  found  to 
be  true  In  all  cases  when  a  guiding  field  was  absent.  Taking 
the  center  of  each  peak  (or  the  twin  peak  In  the  case  of  17a.), 
the  time  of  transit  between  positions  was  determined  for  a 
number  of  runs  and  averaged.  The  resulting  velocities  of  the 
center  of  mass  of  the  plasmold  in  regions  a.  and  b.  were  then 
obtained  by  dividing  11. 5  cm  and  10.5  cm  by  the  respective 
transit  times.   An  upper  bound  to  the  thermal  velocity  was 
obtained  by  choosing  the  point  on  each  curve  which  was  at  half 
maximum  following  the  peak  and  similarly  determining  its  velocity, 
These  velocities  subtracted  from  those  of  the  center  of  mass 
gave  a  measure  of  the  speed  with  which  the  plasmold  was  expand- 
ing.  For  example,  the  center-of-mass  velocities  In  regions  a. 

and  b.  were  v„  =  5<>5  cm/p.sec,  v,  =  1.50  cm/iisec.  The  velocities 
a  D 

of  the  half-maxima  following  the  peak  were  3«8  and  0.95  cm/|isec, 
respectively.  Therefore,  the  expansion  velocities  in  regions  a. 
and  b.  were  at  most  v.   =1.7  cm/p-sec  and  v.^  =  0,57  cm/|xsec, 

A  representative  set  of  data  for  time-of-flight  studies  in 
a  magnetic  field  is  shown  in  Figure  18  (again  signals  are 
negative-going  and  of  arbitrary  amplitude),   Photomultipller 
positions  are  the  same  as  before.  The  double  peak  is  now  not 
only  distinct  in  the  first  position  (Figure  l8a » )  but  it  appears 
to  be  maintained  to  some  extent  during  the  transit  of  the 
plasmold  along  the  tube  (4  and  8  p-sec  in  Figure  l8b.  and  6.5 
and  15.5  usee  in  Figure  l8c.).  This  suggested  that  the  sub- 
bursts  ought  to  be  analysed  separately,  which  was  done.  They 
will  be  identified  as  the  leading  (L)  and  trailing  (T)  peaks. 
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Center-of-mass  velocities  were  calculated  from  assembled 
data  for  both  peaks  L  and  T  in  regions  a.  and  b.  Expansion 
velocity  determinations  for  peak  T  were  readily  performed  as 
In  the  case  of  no  field,  above,  utilizing  the  variation  in 
velocity  of  the  half -maximum  which  followed  peak  T.   In  the 
case  of  peak  L,  the  shape  of  the  curve  did  not  lend  Itself  to 
this  procedure  (e.g,,  using  the  half-maximTom  which  preceeded 
L)  and  there  were  problems  connected  with  relative  rise-times 
of  the  three  photomultipliers.  Therefore,  a  reasonably  repro- 
ducible technique  was  employed  for  estimating  the  change  in  the 
half-width  of  peak  L|  this  result  was  attributed  to  expansion 
velocity  associated  with  temperature. 

Table  1  summarizes  the  results  of  the  photomultiplier 

study,  where  v  ,  v  designate  center-of-mass  velocity  in  regions 
a   b 

a.  and  b,,  and  v,  ,  v.,  designate  expansion  velocity? 

^3    ^  D 


Table 

1 

No  Field 

(L) 

Field 

(T) 

% 

5.5 

5.5 

3o3 

cm/p-sec 

^b 

1.3 

5.3 

1.65 

cm/jisec 

^ta 

1.7 

1,8 

1.9 

cm/M-sec 

^tb 

0.37 

1.9 

0,55 

cm/|jisec 

(V-b)' 

18 

1 

4 

(vta/-tb)' 

21 

1 

12 

The  fifth  and  sixth  lines  represent  the  decrease  in  energy  for 
center-of-mass  velocity  and  expansion  velocity,  respectively. 
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Figure  18 
Photomultiplier  Traces  of  Plasmoid  Transit  with  Field 
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The  following  observations  may  be  made»  First,  without 
field  the  plasmoid  behaves  predominantly  as  a  single  puff  of 
plasma  which  gradually  loses  kinetic  energy  of  both  center-of- 

mass  and  expansion  motion  at  about  the  same  rate.  One  may 

11 
presume  that  residual  gas  within  the  tube  (density  of  1.6  x  10 

-6 

molecules/cc  at  pressul^e  of  5  x  10   mm.  Hg)  barely  contributts 

to  this  deceleration  (mean  twee   path  in  Np  is  about  two  meters). 
However,  the  limitation  of  the  tube  wall  has  become  a  major 
factor  by  the  time  the  plasmoid  has  reached  region  b.   When  a 
guiding  field  is  provided,  however,  the  two  sub-bursts,  which 
appear  to  be  formed  initially,  maintain  their  identity.  The 
behavior  of  L  In  region  a.  with  field  Is  identical  (within 
experimental  error)  with  that  of  the  entire  plasmoid  without 
field  (i.e.,  both  center-of-mass  and  expansion  velocities  are 
the  same).   However,  with  field,  L  is  able  to  maintain  the  mag- 
nitude of  both  velocities  in  region  b.   Conversely,  T  not  only 
starts  (in  region  a.)  with  both  velocities  lower  but  has  lost 
considerable  energy  by  the  time  it  has  reached  region  b. 

Third,  it  is  interesting  to  compare  energy  loss  without  and 
with  field  (the  T  peak)o   As  noted  above,  the  former  case  shows 
equal  energy  loss  for  both  center-of-mass  and  expansion  velo- 
cities, while  in  the  latter  case  energy  loss  associated  with 
the  center  of  mass  is  the  smaller  of  the  two,  and  both  are  less 
than  is  the  case  without  field.  For  one,  a  confining  field 
will  eliminate  energy  losses  to  the  walls.   In  addition,  higher 
density  shortens  the  mean  free  path  for  all  two-body  processes, 
possibly  increasing  the  rate  of  recombination,  and  radiative 
cooling. 
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The  expansion  velocity  of  1.9  cm/M.sec  for  atomic  weight  64 
may  be  associated  with  an  ion  "temperature"  of  120  eV  by 
W:.=  Mv  /2  =  kT.   There  are  several  reasons  for  which  the  true 
kinetic  temperature  n.ay  be  less.   For  one,  the  axial  expansion  of 
the  plasmoid  from  its  center  of  mass  may  be  due  to  a  spatial  varia- 
tion in  initial  axial  velocity  among  the  ions  which  received  the 
initial  impulse.   Second,  the  distribution  of  relative  velocities 
may  be  strongly  anisotropic  at  t  =  0,  and  this  spread  in  axial 
velocities  may  not  have  become  randomized.   Third,  while  it  Is 
convenient  to  assume  that  the  plasmoid  is  macroscopically  neutral, 
this  need  not  be  exactly  true.   Small  deviations  from  neutrality 
would  also  contribute  to  an  increase  in  expansion  velocity,  not 
properly  equatable  to  temperature.   Therefore  w.  =  120  eV  must 
be  considered  an  upper  bound.   This  is  quite  adequate  to  show 
that  the  flow  is  supersonic.   It  i:  of  Interest  to  note,  however, 
that  the  plasmoid  expands  rapidly  to  the  tube  walls  when  no  con- 
fining field  is  applied!  this  may  be  seen  most  clearly  in  Figure 
22b  (vide  infra )  which  must  be  compared  to  Figure  2  for  appropriate 
dimensions.   Here  the  plasmoid  has  a  distance  of  about  one-half 
inch  to  expand  to  the  tube  wall,  and  this  is  accomplished  in  an 
axial  distance  which  is  clearly  smaller  than  two  inches,  suggest- 
ing that  even  radial  expansion  velocity  is  greater  than  one- 
quarter  of  the  axial,  center-of-mass  velocity.   Lack  of  neutrality 
would,  of  course,  be  an  operative  factor  in  this  case,  also. 

The  formation  of  sub-bursts,  L  and  T,  probably  is  due  to 
the  ringing  of  the  plasma  source  capacitor.   If  one  ascribes  the 
velocities  of  5.5  and  3.3  cm/p-sec  to  L  and  T,  respectively,  at 
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the  position  of  the  first  photomultiplier,  8  cm  from  the 
source  mouth,  then  L  required  lo5  ixsec,  and  T,  2 A   p-sec,  to 
travel  that  distance.   Since  Figure  l8a «  shows  L  as  having 
arrived  2.6  |j-s«c  after  (some  arbitrary)  time  zero,  and  T  at 
5.2  tisec,  then  L  started  at  the  source  mouth  at  1.1  ixsec  and 
T,  at  2.8  iisec.   The  two  peaks  were  therefore  separated  a 
time  1.7  M,sec  Initially,  a  frequency  of  6lO  kc,  corresponding 
closely  to  the  measured  620  kc  of  the  capacitor  (Section  IIB). 
3.   Derived  Results;  An  Indication  of  the  quantity 
of  material  released  on  a  single  burst  was  obtained  by  firing 
once  with  the  vacuum  manifold  closed,  and  measuring  the  pres- 
sure.  This  was  more  likely  to  be  a  measure  of  entrapped 
electrode  gas  than  of  metal  ions  but  it  provided  a  rough  figure 
for  test  purposes  I  there  was  no  doubt  that  it  correctly  reflected 
electrode  (plasma  source)  decomposition  since  pressures  obtained 
were  proportional  to  the  square  of  the  plasma  source  voltage. 
A  reading  of  5.8  x  10   mm  Hg  was  obtained  for  a  system 
voltome  of  lj5.8  liters.   The  plasmoid  was  formed  with  a  diameter 
of  about  1.5  cmi  it  passed  the  first  photomultiplier  with  a 
mean  speed  of  4.4  cm/iisec  and  lasted  for  about  12  p.sec  (see 
Figure  l8a.)o   Total  length  was  therefore  about  53  cm  and 
volume,  94  cc.   Ion  pressure  within  this  plasmoid  would  thus 
have  been  (5.8  x  10"  )( 13800/94)  ^  8.5  x  10"^  mm   Hg  at  room 

temperature,  and  the  corresponding  mean  number  density,  n, 

15     / 
utilizing  perfect  gas  law)  was  about  2.7  x  10  ^  ions/ccj  in  the 

17 
94  cc  there  were  2.5  x  10   ions. 

The  method  employed  is  questionable  and  the  result  is  high. 
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A  significant  portion  of  the  plasmoid  may  have  contained  trapped 
neutral  material  which  was  carried  along.   The  result  should  be 
regarded  as  an  upper  bound. 

One  may  test  Its  validity  to  some  extent  as  follows  although 
the  above  objection  would  not  necessarily  be  obviated.   Assiame 
that  the  center-of-mass  velocity  Imparted  to  the  leading  sub-burst 
was  derived  entirely  from  the  magnetic  energy  resulting  from  the 
discharge.   The  Ion  density  at  the  source,  of  this  sub-burst,  will 

not  be  the  same  as  n,  above,  but  higher.   Assume  that  not  more 

17 
than  10   Ions  are  contained  In  this  sub-burst.   Cross-section 

will  be  that  of  the  center  electrode  (location  of  highest  current 
density)  where  d  =  0.33  cm.   A  sub-burst  lasts  about  1.7  iisec 
(vide  supra ) ;  at  v  =  5.5  cm/M.sec  this  amounts  to  9-5  cm  of  length. 
Volume  Is  therefore  0.8  cc  and  density  Is  1.2  x  10  '  lons/cc, 
at  the  source.   The  kinetic  energy  density  may  be  written  as 
pv  /2  where  p  is  the  mass  density  and  v  is  the  resulting  velocity. 
This  may  be  equated  to  the  magnetic  energy  density  p.  H  /2.   The 
value  of  H  is  not  known  directly  for  the  plasmoid  discharge,  but 
one  may  substitute  l/i  where  i  is  some  characteristic  distance. 
Since  the  center  electrode  diameter  is  about  0.33  cm,  t   will  be 
taken  as  0.33'ir  x  10   m.   A  rough  measure  of  I  (at  maximvun  dis- 
charge) is  given  by  (3/2)(q/At)  where  q  is  the  initial  charge  on 

the  capacitor  and  At  is  the  time  for  a  quarter  wave  (in  this  case, 

7  k 

4.1  X  10  '  sec).   At  12  kV  and  0.55  ^iF,  I  =  2.4  x  10  amp.   For 

17  — P      ^ 

n  =  1.2  X  10   lons/cc,  p  =  1.3  x  10~  kg/m   (taking  a  mean  atomic 

weight  of  64).   Therefore,  one  obtains 
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I  2A   X  10^ 
\  0.53Tr  X  10"^' 

s  2ol  X  10  m/sec 
-  2„1  cm/p.sec, 

which  is  to  be  compared  to  the  measured  velocity  of  5«5  cm/^isec| 
the  2o 5-fold  discrepancy  reflects  the  uncertainlty  in  the 
computation. 

These  numbers  may  be  used  for  another  purpose,  to  estimate 
an  upper  bound  for  the  self-field  created  by  any  possible  plas- 
moid  current So  The  simplest  assumption  would  be  that  the 
initial  discharge  current  and  consequent  field  is  larger  than 
any  subsequently  developed  in  the  plasmoid.  The  induction 
created  in  this  case  is  then  given  directly  by  B  =^  ^q^'^^  ~ 
3  weber/m  ,   However,  the  plasmoid  cro§s-section  increases 
rapidly  by  a  factor  of  (lo5/0,33)  ^21,  so  that  the  plasmoid 's 
effective  field  would  be  well  below  that  of  the  solenoid. 

It  is  of  interest  now  to  consider  the  nature  of  the  plas- 
moid discharge  (with  particular  attention  to  ion  and  electron 
temperatures  and  degree  of  ionization)  during  the  time  that 
current  is  flowing  between  the  electrodes  and  shortly  thereafter, 
a  period  of  the  order  of  4  M,sec.   Continuous  arcs  have  been 
studied  extensively  as  part  of  the  general  subject  of  gas 
discharge,  with  Edels"^^  making  a  particular  subject  of  the  study 
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of  the  temperature  of  irdivldual  species  within  the  discharge. 
The  problem  here  is  a  rather  different  one,  however,  in  that 
times  of  interest  are  extremely  short  and  the  discharge  is  formed 
along  a  surface  in  initially  high  vacuxom  conditions.   Transient 
vacu\am  sparks  have  been  described  by  Boyle,  et^  al^=   For  small 

(order  of  mm)  gaps  at  high  voltage  (>  10  )  and  the  problem  has 

12 
been  examined  with  respect  to  arcing  in  switches  by  Suits 

some  years  ago.   Here  the  discharge  is  probably  initiated 

largely  by  electron  current,  with  initial  energies  possibly 

comparable  to  the  gap  voltage,  and  subsequent  vaporization  of 

the  electrodes  from  which  a  low  degree  of  ionization  results. 

With  the  formation  of  sufficient  ionized  vapor  a  "spark  channel" 

is  created   through  which  high  currents  may  now  passj  these 

processes  can  occur  in  times  substantially  shorter  than  1  iisec. 

There  is  little  presently  known  about  the  relative  tempera- 
tures of  ions  and  electrons  during  this  early  period.   Presumably 
the  electron  temperature  is  initially  very  high|  in  colvunn 
discharges  Massey   reports  w  »  w„ ,  as  does  Edels   for  the 
particular  case  of  low  pressure  discharges,  but  the  times 
involved  are  very  much  longer  than  the  duration  of  the  experi- 
ments reported  here. 

It  seems  generally  agreed  that  one  can  have  some  form  of 
temperature  equilibration  in  the  order  of  milliseconds j  distant 
coulomb  collisions  are  effective  for  this  purpose  with  high 

initial  electron  temperatures,  as  may  be  seen,  for  example,  by 

IS 
the  equipartition  times  given  by  Spitzer,  ■"     He  gives  the 

relation 
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n,z  z- 


as  the  time  for  equlpartltion  of  energy  for  a  test  particle  of 
mass  number  A  among  field  particles  A,,   In  (1),  A  is  the  ratio 
of  cut-off  distance  for  interactions  (Debye  shielding  dis- 
tance) to  the  distance  of  closest  approach,  and  InA    is  of  the 
order  of  tenj  T  is  in  degrees  Kelvin.   In  this  case,  one 
wants  the  equipartition  time  for  high  energy  electrons  found 
in  the  initial  discharge  to  reduce  to  the  ion  temperatures 
(field  particles).  Thus 

c-  A.A   /w  \  3/2 

since  w  /A  »  w^j/A.  if  any  reasonable  value  of  w  is  chosen 
(w  in  eV).   If,  for  example,  one  chooses  w  =  120  eV,  and 
w  =  1200  eV,  n^  =  1.2  x  10  ,  and  assume  a  Cu  plasma  (A.  =  64 )j 
substitution  in  (2)  yields  times  of  the  order  of  500  ixsec  for 
z.  =  1  or  50  M-sec  for  z.  =  4  = 

Ionization  times,  however,  are  considerably  faster  (as 
will  be  shown  shortly)  and  one  may  suppose  that  the  speedy 

attainment  of  low  electron  temperatures  noted  by  Meek  and 

1"5 
Craggs  ^   and  others  is  due  to  a  rapid  increase  in  ionization 

s 
level,  with  the  initial  hot  electrons  losing  their  energy  by 

ionizing  and  heating  bound  electrons j  the  ion  temperatures 

would,  of  course,  be  largely  unaffected  by  this  process. 

Suppose,  for  example,  that  during  the  early  stage  of  the  arc 
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formation  w^  had  reached  a  level  of  several  kilovolts  while 
e 

the  copper  vapor  was  still  largely  singly  ionized,  and  w.  = 
120  eV  (these  conditions  will  be  taken,  for  the  moment,  as 
entirely  hypothetical  in  order  to  serve  the  computations  which 
follow).  The  free  electrons  present  would  proceed  to  increase 
the  ionization  level,  dropping  in  temperature,  while  losses 
through  coulomb  interactions  with  the  ions  would  be  negligible. 
The  question  arises,  where  would  an  equilibri\im  of  some 
variety  be  achieved? 

These  plasmoids  are  optically  thin  so  that  the  Saha  equa- 
tion  is  undoubtedly  inapplicable.   In  1952,  Elwert   developed 
an  equation  for  collision  ionization  -  photorecombinatlon 
equilibrium  in  the  sun's  corona: 

!l±l   -  fl2  -  8  J  X  105  '2  ''''  ''hV  '''p'-V^'     (3, 

"i   "  «2a  -  S-'  ^  ^°  17  ?^  [III      (Vj    .  (5) 

where 

Ni+1 

'  'u  ~   relative  prevalence  of  ion-(l  +  1  ionized)  vs 

"i 

ion-(i  ionized) , 

3^2  =  electron  collision  probability  coefficient, 

Q22  =  photorecombinatlon  coefficient, 

f-j^/fp  =  uncertainty  factor  =  9/8, 

L  =  number  of  electrons  in  outemnost  occupied  shell, 

n*  =  principal  quantvim  number  of  outennost  shell, 

g  =  recombination  probability  factor  for  n*  or  higher  shells, 

a  function  of  I„(n*)  /w  given  in  graphic  form  by  Elwert. 
n        e 

If  one  takes 
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where  Z  -.-.  is  the  effective  nuclear  charge  seen  by  the  i 
electron  and  I„  is  the  ground  state  ionization  potential  of 
hydrogen  (13.53  eV),  and 

I. 
s=rA,  (5) 


then  combining  (3),  (^),  and  (5)  yields 


^  .  9.4  X  105  ^  ^^f  s-5  exp(-s).  (6) 

Suppose,  now,  one  chooses  w  =  120  eV  at  equilibrivun  (as 
was  just  hypothesized)  and  enquires  as  to  the  most  prevalent 
ionized  species  of  copper  (Z  =  29).  This  may  be  determined 
by  a  trlal-and-error  process.   For  example,  one  finds  that 
for  i  =  z  =  15 

Cu(XV)2   (ls)^(2s)^(2p)^(3x)^ 

N,^/N, ^  =1.4  (the  sub-states  of  the  n*  =  3  electrons  are 
lb   15 

uncertain  here).   In  this  calculation,  take  Z^^f.  =  17-5,  n*  =  J>, 
L  -  4,  from  which  s  =  3.85,  s(n*)^  =  34.6,  g  =  4.3  (from 
Elwert's  graph).   For  1  =  z  =  I6 

Cu(XVl);   (ls)^(2s)^(2p)^(3x)^ 
^17^\e   "  °-5  (^eff  ^  ^^°5,  n*  ==  5,  L^  =  3,  from  which  s  =  4.3, 

-  61  - 


s(n*)^  =  38.5,  g  =  ^.5).   Therefore,  at  w  =  120  eV  the  most 
prevalent  specie  Is  Cu(XVI),  i.e.,  z  =  I6. 

Suppose  the  Initial  conditions  before  equilibration  were  as 
stated  above,  namely  w  >  w. ,  Cu(l)  present.   Then  each  electron 
would  have  to  produce  I5  ionizations,  a  total  energy  of 
(approximately) 


16    2 


This  may  be  estimated  as  2100  eV.   In  addition,  I90O  eV  must  be 
distributed  among  I6  electrons  to  produce  120  eV  each.  The 
initial  electron  must  then  have  had  about  4  keV.   This  is 
admittedly  high. 

It  is  necessary  also  to  show  that  the  relaxation  time  for 

this  process  is  sufficiently  short,  and  this  may  be  done  utiliz- 

17 
Ing  an  expression  developed  by  Knorr  '  for  corona  equilibrium. 

Using  symbols  as  defined  above, 

6.5  X  lO^Csw^)^-^^  exp(s) 

*  = ri — ^ (^) 

n  e 

where  n  ,  the  electron  number  density,  is  given  by  zn  .   For 
"  i 

Cu(XVl),  s  =  4.5,  w^  =  120  eV,  L^  =  3,  and  n^  =  (l6)(2.7  x  lO""-^). 
Substitution  of  these  values  in  (7)  yields  t  =  4.3  |j.sec  which 
is  of  the  order  of  the  experiment  time.   It  is  likely  to  be  an 
upper  bound,  however,  since  (?)  is  constructed  primarily  for  an 
examination  of  equilibration  time  for  electron  energies  close  to 
final  equilibritmi.  But  in  the  process  suggested  here  the  initial 
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electron  temperature  would  be  considerably  higher,  and  equilib- 
rium would  be  attained  at  a  faster  rate. 

An  initial  electron  temperature  of  4  keV  is  difficult  to 
accept,  although  there  might  well  be  a  small  proportion  of 
"runaway"  electrons  produced  in  the  initial  discharge.   However, 
their  thermalization  would  depend  upon  a  rapid  increase  in 
electron  scattering  as  the  spark  becomes  pinched  prior  to  its 
propulsion  down  the  tube.   It  is  likely  that  the  electron 
temperature  is  initially  one  to  two  orders  of  magnitude  smaller. 
Consider,  for  example,  a  corona  equllibrlxim  at  w  =12  eV.  For 
1  =  z  =  6 

Cu(VI) s   (ls)^{2s)2(2p)^(3s)^(3s)^(3d,4s)5 

N^/Ng  =5.4  using  (6).   The  quantum  states  of  the  five  outer 
electrons  are  uncertain  but  the  5d  and  4s  levels  should  be 
sufficiently  close  so  that  ionization  from  either  is  equivalently 
weighted.  For  this  computation,  take  Z  »^  =  7,5,  n*  =  5.3, 
Ljj  =  5,  from  which  s  =  5»75,  s(n*)^  =  6:5,  g  =  5»3.   For  i  =  2  =  7 

Cu(VIl)s   (ls)2(2s)2{2p)^(5s)2(3p)^(3d,4s)^ 

Ng/N^  =  0.9  (taking  Z^^^  =  8.5,  n*  =  3^5,  L^  =  4,  from  which 

s  =  6.75,  s(n*)  =  81,  g  =  5.7)-  Thus,  the  prevalent  species  at 

w^  =  12  eV  is  Cu(VIl),  i.e.,  z  =  7- 

If  initial  conditions  before  this  equilibrium  were,  againp 
that  w^  >  12  eV,  Cu(l)  prominant,  then  each  electron  would  have 
to  produce  6  ionizations  for  a  total  energy  of  (approximately) 
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"  z»2  (n*)^ 


This  may  be  estimated  as  175  eV.   In  addition,  85  eV  must  be 
distributed  among  7  electrons  to  produce  12  eV  each.   The  initial 
electron  must  then  have  had  about  260  eV,  or  about  twice  the  ion 
temperatxire.   This  is,  perhaps,  a  more  reasonable  initial 
condition. 

The  relaxation  time  for  this  process  should  be  tested, 
also.   Using  the  above  parameters  for  Cu(VIl),  with  n  = 
(7)(2.7  X  lO'"'^),  (7)  gives  t  =  5.^  usee.  This  is  again  of  the 
time  order  of  the  experiment  and  should  be  considered  an  upper 
bouad  for  the  reasons  previously  stated. 

For  the  purposes  of  understanding  the  origin  of  the  high 
ion  temperatures  suggested  by  the  data,  one  may  speculate  on 
a  mechanism  whereby  the  plasmoid  might  acquire  a  thermal  energy 
comparable  to  its  energy  of  axial  motion.   Experience  with 
discharge  phenomena  suggests  that  the  initial  ion  temperature 
may  have  been  no  more  than  a  few  electron  volts.    However, 
it  is  probable  that  the  spark  channel  pinches  during  current 
passage.   The  effect  upon  the  sub-burst  is  to  blow  it  outward, 
down  the  tube,  and  might  be  likened  to  the  sub-burst  being  hit 
Tf>y  a  fa  St -moving  wall.   This  creates  a  shock  wave  in  the  gas, 
ind  a  significant  fraction  of  the  axial  kinetic  energy  acquired 
from  the  pinch  can  thus  be  thermallzed.   In  the  limit  of  a 

strong  shock  in  a  polytropic  gas  this  fraction  can  be  approxi- 

18 
mated  theoretically.    One  writes 
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where  M  is  the  initial  Mach  niunber,  the  subscripts  refer  to  the 
states  on  either  side  of  the  shock,  and 


2   7-1 

^      Y  +  1 


For  this  discussion  take  y  -   5/3,  M-  =  1/2.   If  one  assumes  (as 
has  been  done  here)  that  the  axial  velocity  imparted  to  the 
plasmoid  is  approximately  that  of  the  moving  wall  (the  shock), 
then  M  =  v  /c,  with  c  the  sound  speed  in  the  unperturbed  spark 

O       a 

channel.   But  one  has 


c^  =  yR'T^  , 


where  R'  is  the  gas  constant  per  unit  mass,  which  results  in 

2 

p         P   V 

_  ^Jt  (1  +  t^  )  a 
^1  "      7     R^~  ° 


Now  (3/2)R'T,  is  the  final  thermal  energy  per  unit  mass, 
acquired  from  interaction  with  the  shock,  of  which  the  ions 
acquire  half,  i.e.,  T^  =  (2/3)v^g/R'.   Therefore,  one  obtains 


V^2  ^  3^2(1  +  u^)  _  0  28. 
2t 


Insofar  as  the  endothermic  reactions  behind  the  shock  have 
been  neglected  (increase  in  degree  of  ionization)  this  is  an 
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overestimate.   Experimentally,  for  the  leading  sub-burst 
(v^/v^)2  =  0.11. 

One  may  also  show  that  bremstrahlung  is  not  an  appre- 
ciable energy  loss  for  these  plasmoids.   Taking  the  power 

IS    /    \ 
loss  per  unit  voliome  given  by  Spitzer  -^  as  (emu) 

1.4  X  10-2-^  z\2  T^l/2  , 

one  obtains  for  the  energy  loss  per  electron  during  the 
course  of  the  experiment  (12  usee) 

Aw  =  1.05  X  10"^-^  z^n^  T^^^'^   eV. 

15 
An  upper  bound  is  obtained  by  taking  z  =  16,  n.  =  2.7  x  10  ■^, 

T  =  1.2  X  10  ,  from  which  Aw  =  0.8  eV  per  electron. 

With  these  possible  values  of  z  (depending  upon  assximed 
initial  conditions),  theoretical  estimates  may  now  be  made  of 
the  plasma  drift  velocity,  and  compared  to  the  observed  expan- 
sion of  the  plasmoid  flow  in  region  c,  as  described  in 
Section  1,  above.   There,  a  one  percent  lateral  drift  was  noted. 
If  V|^  =  5.5  X  10  cm/sec,  then  the  expectation  is  that 
v^  =  5  X  10  cm/sec. 

There  are,  unfortunately,  two  non-equivalent  drift  proposals 
from  which  to  choose  at  the  present  time.   The  now  classical 
interpretation  is  that  complete  transverse  confinement  is  depen- 
dent upon  infinite  plasma  conductivity,  with  lateral  drift  being 
proportional  to  E  x  B  (the  electric  field  resulting  from 
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resistance  within  the  plasma),  the  result  of  electron-ion 
collisions.   This  leads  to  an  expression"^^  (emu) 

Conversely,  Bohm,  Burhop,  and  Massey^  have  proposed  (without 
proof)  that  drift  in  the  presence  of  a  magnetic  field  is  a 
"drain"  mechanism,  a  consequence  of  the  crossed  electric  fields 
caused  by  plasma  oscillations,  and  that  the  drift  velocity  is 
given  by  ^   (emu) 

V,  .M^-io^l^l  .  (9) 

e 

It  is  interesting  to  note  that^  aside  from  a  difference  in 
dependence  on  field  strength,  (8)  is  proportional  to  z/t  ' 

while  (9)  depends  implicitly  on  T/z  (noting  that  n  =  zn. ) . 

Purported  experimental  verifications  for  both  (8)  and  (9)  have 

20 
been  recently  presented,  respectively,  by  Simon   and  by 

Lehnert . 

Proceeding  now  to  apply  these  expressions  to  the  present 

case,  an  estimate  of  V9   is  required.   An  upper  bound  for  the 

average  pressure  is 

p  =  n^kTj^  . 

If  one  assumes  a  circular  cross=>section  for  the  plasmold  and  a 
linear  (radial)  pressure  gradient,  then  it  can  be  easily  shown 
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that 


Pmax  =  5P  =  5"l^i  • 


The  plasmoid  diameter  is  about  I.5  cm,  from  which 

\Vp\    -^-   S^T^  .  (10) 

Consider,  first,  the  electron-ion  collision  mechanism. 
With  ln^    =  10,  B  =  1.1  x  10  ,  n^  =  2.7  x  lO'''^,  T^  =  1.2  x  10^, 
one  obtains  from  (8)  and  (10)  v  =  I.5  x  lO^z  cm/sec.   If  a 


very  high  initial  electron  temperature  is  assiomed,  as  above, 

4 


perhaps  z  =  I6,  and  v  =  2.4  x  10  cm/sec.   This  is  not 


unreasonable  in  comparison  to  an  experimental  value  of  5  x  10 
cm/sec.  But  it  is  seen  that,  for  this  mechanism,  any  ass\imption 
of  lower  z  increases  the  discrepancy. 

Employing  the  Bohm  proposal,  a  combination  of  (9)  and  (10), 

with  parameters  as  before,  yields  v^^  =  (2.4  x  lOvz)  cm/sec; 

^  4    / 

If  z  =  16,  Vj_  =  1,5  X  10  cm/sec, 

ifz=7,  v^=5.5xlO  cm/sec, 
the  latter  being  closer  to  observation.   It  is  seen,  incidentally, 
that  the  Bohm  drift  velocity  will  not,  in  general,  be  much  the 
larger  in  strong  magnetic  fields  except  in  the  special  case  of 
fusion  interest,  namely,  a  hydrogenic  plasma  at  moderate  density 
and  very  high  temperature. 

One  may  now  calculate,  also,  a  range  of  possible  ion  Larmor 
radii.   Using  (emu) 

'^^ta 
^   Bze  ' 
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a  measure  of  r  may  be  obtained  for  the  a.  and  c,  regions.   If 
M  =  64  a.m.u.  in  grams,  v.^  =  1,9  x  10  cm/sec,  B  =  1.1  x  10 
emu,  z  =  7,  then  r  =  0,17  cm.  This  will  be  reduced  to  r  = 
0.08  cm  for  region  b.   Any  higher  value  for  z  and/or  lower 
ion  temperature  will  reduce  these  figures  accordingly.   In  any 
case,  they  are  substantially  smaller  than  the  plasmoid  dimen- 
sions and  any  distance  of  appreciable  variation  in  Bo 

Another  question  of  interest  concerns  the  interaction 
between  the  plasmoid  ions  and  the  "mirror"  separating  regions 
a.  and  b.  in  the  normal  Differential  Solenoid,  The  effect  of 
a  magnetic  mirror  geometry  is  well  known  in  problems  of  plasma 
confinement;  it  depends  upon  the  fact  that  the  magnetic  moment 
of  the  individual  ion  in  a  helical  path,  |i  -  eVj^r/2,  is  an 

adlabatic  invariant  of  the  motion  (if  B  is  constant  over  r), 

z 

so  that  a  positive  B  gradient  behaves  like  a  potential 
barrier.   The  particle  thus  gains  transverse  kinetic  energy, 
w  ,  as  the  Larmor  radius  decreases,  which  by  conservation  it 

obtains  from  the  axial  kinetic  energy.   If  the  mirror  ratio, 

22 
R  =  B^/b  ,  is  sufficiently  large,  the  particle  is  reflected j 

the  critical  angle  for  reflectlon'^^  is  sin^  0  =  l/R.  This 
analysis  can  be  extended  to  a  dynamical  plasma,  i.e.,  a  plas- 
moid with  a  cm.  velocity  along  the  axis,  with  rather  interesting 
results.   In  particular,  in  the  limit  of  zero  cm.  velocity  it 
will  be  shown  that  the  reflection  coefficient  is  F^  =  (l-l/R)-"-/^. 
This  total  reflection  coefficient  is  the  fraction  of  the  entire 
ion  population  in  the  plasmoid  ultimately  reflected,  and  should 
not  be  confused  with  the  1=1/R  of  Spltzer,  which  is  the 
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fraction  of  the  Incident  ions  reflected  per  unit  time. 
To  reformulate,  total  kinetic  energy  is  conserved; 

dWj,  +  dWj_  =  0,  (11) 

where  w,,  and  w,  refer  to  axial  and  rotational  kinetic  energy. 
If  \i   is  an  invariant,  then  one  may  write 

dWj^  =  \x6B.  (12) 

Combining  (11)  and  (12),  one  obtains 

where  w,j  and  B  may  be  assumed  (for  simplicity)  to  have  only  z 
components  which,  in  turn,  have  only  z  gradients.   The  first 
term  of  (13)  may  be  integrated  as  a  force  decreasing  the  axial 
velocity,  so  that 

dv    dw 


dt    dz 


•.  § .  w 


In  the  simplest  case,  all  differentials  can  be  taken  as  linear 
variations,  yielding 

Av.        .T3 
"SF^'-^'i-  (15) 

Take  the  plasmoid  to  consist  of  N  Ions  (and  associated 
electrons)  having  an  isotropic  Maxwelllan  velocity  distribution, 
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f(v),  with  respect  to  the  cm.,  and  with  a  cm.  speed  Zf, 
the  direction  of  which  defines  the  axis.   A  fraction  of 
the  ions,  dN/N,  will  possess  a  speed  in  the  range  v  and 
v+dv,  and  thus  will  lie  on  the  surface  of  a  sphere  of  radius 

V  in  velocity  space,  with  axis  parallel  to  the  cm.  motion. 
In  the  moving  frame  of  the  cm.  motion,  furthermore,  the 

V  of  every  ion  contained  in  dN  will  lie  on  some  cone  of 
revolution  with  angle  a  to  the  axis  (O  <  a  <  tt).   We  now 
enquire  what  fraction,  f_^,  of  the  dN  ions  with  velocity 

V  will  be  reflected  as  a  consequence  of  their  orientation, 
when  the  plaamoid  enters  the  magnetic  mirror. 

First,  take 

v^  =  V  sin  a,  (l6) 

V,,  =  V  +  V  cos  a.  (17) 

The  expression  for  |jl  can  now  be  written 


M-  = 


^  sin  a  (18) 


from  (16).   Combining  (15)  with  (I8),  one  obtains 


2MAv,,Az 
^^^  <=^-   -   evrAtAB  *  ^^^^ 


If  the  particle  is  to  be  stopped,  Av„  =  -v„ .   If  the  particle 
is  experiencing  uniform  deceleration,  Az/At  =  V||/2.   Also,  one 
may  write  AB  =  (R-l)B^.   Substituting  in  (19),  one  obtains 

9 
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M(v„)^  (20) 

siJ^  a  =  evB  (H-i)  • 

a 

However,  there  is  also  the  relationship 

(21) 
B  er  =  Mv,  =  Mv  sin  a.  ^ 

a       -^ 

Combining   (1?),    (20),   and   (21)   yields 

V/v  =    (R-l)-*"^^   sin  a   -   cos  a.  (^^^ 

(The  result  may  he  obtained  also  by  transforming  sin^  Q   =  l/R 
to  a  coordinate  system  moving  with  the  cm.)   This  result,  (22), 
may  be  plotted  as  the  speed  ratio,  S  =  Vu,  versus  a,  a  family  of 
curves  in  R.   Figure  19  shows  curves  for  several  R,  as  noted.   It 
is  convenient  to  solve  (22)  for  cos  a,  yielding 

r     2     n    U/2 
COS  a  =  -  i^  ±  f l|  (^)  -lltH-llj 


1/2  /^,N 


where 


i         .  -  n  .  1  .  (24) 

R 


a  H^  ,     ^  -  1  -  R 


This  double-valued  function,  (25),  defines  the  angular 
limits  ai  and  a^,  between  which  an  ion  with  velocity  v   must  lie 
in  order  to  be  reflected;  if  one  takes  the  positive  sign  to 
correspond  to  a,,  then  0  <  a,  <  a^  <  -   For  a  given  mirror  ratio, 
R,  and  axial  speed,  U,  a  pair  of  angles  for  reflection  are 
defined  only  down  to  some  minimum  value  of  v,  at  and  below  which 
no  orientation  will  provide  reflection.   One  sees  from  (25)  that 
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l/p 
this  mlnlmxim  value  is  v  =  a  '^   for  which  a  single  a  is  defined, 

cos  a  =  -  1/R   .   Every  spherical  velocity  shell  with  v  >  a  '^ 

will  contain  a  band  of  ions,  bounded  by  latitudes  at  a-,  and  a^, 

which  will  be  reflected.   The  fraction  which  these  reflected 

ions  bear  to  the  total  ions  of  any  given  shell,  dN,  is  Just  the 

ratio  of  the  bounded  area  on  the  shell  to  the  total  area.   One 

may  show  readily  that  this  ratio  is  (cos  a,  -  cos  ap)/2<, 

Therefore, 


f/  =  f  1  =  0  I   V  <  a^^2       (25) 


=  f^  =  (cos  a^   -   cos  a2)/2 

1/2 
p)] 

V 


[b(l  -  H_)]  .  V  >  a^/2       ^26) 


based  upon  the  above  arguments  and  (23). 

The  Maxwellian  velocity  distribution  is 


3/2 
^  =  ^ir  (^)    exp  (-  MvV2kT)v2dv 


where 


^1/2         ^     ^ 
=  (l6cVTr)    exp  (-cv^)v^dv,  (27) 


c  =  M/2kTo  (28) 


The  total  fraction,  F  ,  of  the  ions  reflected  is  then 


"r^ 


r  f „^dN 
r 


N 


N 


■I 


fgdN 


1/2     N   * 
N(v  >  a^^'^) 


-   7^ 


(29) 


from  (26).   Combining  the  integral  of  (29)  with  (27)  yields 


r  J 


rv>n   a  n^'^^fi^c^//^   ,   2v  2^ 
lD(l =■)]   [— - — J    exp(-cv  )v  dv 

1/2      v"^       ^ 


*bc511/2  ^ 


— / 


(v^  -  aJ^/2  exp(-cv^)d(v^)      (^0-) 


1/2 
=  b  ^   exp  (-c/a)o  (30) 


(By  adopting  the  transformation  y  =  (v^  -  a)"^^^,  (30')  may  be 
converted  to  integrable  form.)  Now,  applying  (24)  and  (28)  to 
(30)  yields 

Fp  =  (1-1/R)^/^  exp  (-Ml?/2kTR),  (31) 

Now  apply  this  result  to  the  present  data.   A  lower  bound 
for  the  ratio  of  axial  to  thermal  energy  (vide  supra )  for  the 
leading  peak,  L,  is  given  by  (Vg/v^^)^  =  9.3,  and  this  may  be 
equated  to  (2/3) (Maf/2kT) .   The  mirror  ratio  was  R  =  2.5/I.I  =2.3- 
Therefore,  from  (31),  F^  =  O.OOI7,  or  essentially  no  portion  of 
L  was  reflected.   A  similar  analysis  for  the  trailing  peak,  T, 
gives  F^  =  0.11,  which  is  substantially  larger,  but  it  is  prob- 
lematical whether  the  effect  was  photographically  recordable. 
If  one  associates  with  the  T  peak  the  hazy  light  surrounding  the 
central  beam  in  region  a.  of  Figure  l4c.,  then  it  does  appear  as 
if  a  noticeable  fraction  of  these  ions  is  missing  in  region  b. 

There  is  a  possible  extension  of  these  calculations  which 
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may  be  suggested  here.  The  above  analysis  was  designed  to 
investigate  the  possibility  of  stopping  completely  certain 
ions  with  appropriate  velocity  orientation.   In  fact,  even  If 
an  ion  is  not  stopped  it  will  experience  some  deceleration  by 
the  mirror.   A  more  general  solution  might  thus  be  made  to 
yield  the  redistribution  of  velocities  occasioned  by  mirror 
transit  of  the  plasmoid,  based  on  an  initial  isotropic 
Maxwelllan  distribution. 

C.   Flasmold  Flow  In  the  Shorted  Differential  Solenoid; 
The  axial  values  of  magnetic  induction  are  given,  in  each  case 
below,  for  regions  a,  and  c.  and  for  region  b,  (see  Figure  7). 
The  field  Intensity  in  region  b.  decreased  radially  about  25 
per  cent  at  the  wall,  as  is  described  in  Section  IIE2^  while 
the  Jumpers  contributed  a  relatively  strong  transverse  field. 
The  source  was  fired  at  12  kV,  producing  40  joules o  The  pic- 
tures which  follow  are  again  individual  bursts,  taken  with 
open  shutter  and  appropriate  apertures,  and  are  representative 
and  reproducible.   Synchronization  was  arranged  to  place  the 
plasmoid  burst  approximately  10  cm  into  region  a.  at  maximum 
field. 

Figure  20  shows  the  flow  pattern  of  the  plasmoid  in  two 
camera  views.   In  20a,,  is  seen  the  position  of  the  jumper,  with 
the  wires  of  the  second  jumper  cut  and  extending  from  behind  the 
solenoid.   There  is  clear  colllmation  into  region  b.,  where  the 
plasmoid  immediately  expands  to  the  walls  of  the  tube,  as  do 
the  field  lines,  and  is  then  recollimated  in  region  c.  (Figure  20b.) 


-  76  - 


Considerable  turbulence  is  sxiggested  by  Figure  20c,,  with  a 
noticeable  dark  area  along  the  bottom  of  region  b.  This  is  the 
side  on  which  the  Jumper  is  located i  B  at  the  wall  in  the  center 
of  region  b.  is  approximately  0.3  weber/m  while  the  transverse 
field  from  the  jvimper  is  about  0<,5  weber/m  ,  so  that  B  has  a 
strong  azimuthal  orientation. 

To  investigate  this  effect  further,  a  second  Jumper  was 
placed  across  the  opposite  side  of  region  b.  This  would  be 
expected  to  halve  the  transverse  field  toward  the  edge  of  the 
tube  near  the  first  Jumper  and  spread  it  more  uniformly  across 
region  b.  The  somewhat  surprising  result  is  shown  in  Figure  21, 
where  the  incoming  plasmoid  has  been  diverted  to  one  side  of 
region  b.   However,  it  will  be  noted  that  the  emergent  beam  has 
remained  unaltered. 

D.   Plasmoid  Flow  in  the  Curved  Solenoid;  This  geometry  was 
shown  in  Figure  9,  Section  IIF,  and  discussed  there o   The  mean 
value  of  the  magnetic  induction,  B  ,  is  given  in  each  case  below, 
as  determined  by  the  initial  bank  voltage «   The  plasmoid  source 
was  fired  at  3,  6,  and  12  kV,  yielding  2.5,  10,  and  4o  Joules, 
respectively,  as  noted.   Photographs  were  taken  with  shutter  open, 
at  appropriate  aperture,  of  single  bursts,  except  as  noted.   The 
results  shown  here  are  representative  and  reproducible » 

Figure  22  shows  plasmoid  flow  for  two  values  of  mean  B_  and 
plasmoid  at  40  Joules.  The  interesting  result  is  seen  that 
doubling  of  field  intensity  appears  to  cause  no  alteration  in  the 
flow  pathi  in  both  cases  the  plasmoid  is  initially  collimated  but 
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Figure  20 

Plasmoid  Flow  in  the  Shorted  Differential  Solenoid. 
I.  Single  Jumper 
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Figure  22 

Plasmoid  Flow  in  the  Curved  Solenoid, 
I.  Plasmoid  at  L|.0  joules 
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moves  Into  the  tube  wall  at  about  the  sixth  (visible)  solenoid 
turn,  evidencing  little  transverse  confinement.  Beyond  this 
point,  the  plasmoid  presumably  flows  along  the  outer  inside  wall 
of  the  tube.  The  fact  that  doubling  the  field  intensity  pro- 
duces no  visible  effect  suggests  that  the  overriding  processes 
contributing  to  non-confinement  are  not  single  particle  but 
collective,  i.e.,  the  density  is  sufficiently  high  so  that 
single  particle  mechanisms  are  not  substantially  operative. 

With  this  in  mind,  the  experiment  was  repeated  with  the 
plasma  source  capacitor  energized  to  half  the  voltage,  6  kV, 
producing  a  10  Joule  burst.  This  was  expected  to  produce  some 
unspecified  decrease  in  center-of-mass  velocity  and  temperature- 
in  addition  to  quartering  the  density.  The  results  are  seen  i* 
Figure  23.   One  now  notes  a  distinct  Improvement  in  transverse 
confinement  at  1.4  weber/m  and  a  further  improvement  when  the 
field  was  doubled. 

Finally,  the  source  capacitor  voltage  was  again  halved, 
with  the  results  shown  in  Figure  24.  Here  confinement  is  quite 
good,  although  a  portion  of  the  residual  plasma  continues  to 
travel  (and  recombine)  along  the  outer  wall. 

Perhaps  the  least  sophisticated  analysis  that  could  be 
applied  to  these  results  would  be  that  of  comparing  field  and 
plasmoid  energy  density  since  confinement  implies  that  the 
former  exceeds  the  latter,  i.e.,  p  <  1.   For  the  case  of  the 
plasmoid  fired  from  the  source  energized  to  40  Joules  (the  usual 
arrangement  in  this  study)  mass  density  has  been  estimated  at 
2.9  X  10~  kg/m-^  and  maxlmiom  velocity  at  5.5  x  10  m/sec,  from 
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which  the  plasmoid  translational  energy  density  is  about 
4.4  X  10  Joules/m  .   With  magnetic  induction  at  1.4  weber/m  , 
field  energy  density  is  somewhat  higher,  at  7-8  x  10-^  joule s/m^, 
yet  confinement  was  not  observed.   In  addition,  some  visible 

improvement  should  have  resulted  from  increasing  field  energy 

ft  "5 

density  to  3.1  x  10  joules/m  (doubling  solenoid  voltage), 

but  Figures  22c.  and  d.  suggest  otherwise. 

At  the  other  extreme  one  knows  theoretically  that  single 

particles  should  experience  transverse  confinement  under 

appropriate  conditions,  namely  those  which  result  in  drift 

velocities  substantially  smaller  than  the  axial  velocity. 

Single-particle  drift  in  this  geometry  will  result  both  from 

a  transverse  field  gradient,  V|_B,  (such  as  is  present  in  a 

curved  field)  and  from  centrifugal  force.   These  two  velocities 

are  given,  respectively,  by 

V.'  =  -^  (-^  =  -)  (32) 

M    2R       ^  B    R^  ^-^    ' 


and  by 


^a"  =  w-  '  ('') 


where  r  is  the  Larmor  radius,  R  the  radius  of  curvature  of  the 
solenoid,  oo  the  cyclotron  frequency,  and  v_^  and  v,,  the  transverse 
and  longitudinal  ion  velocities,  respectively.   Of  these  two, 
V,"  is  much  the  larger  for  the  present  experimental  parameters 
(and  v^"  for  electrons  is  negligible).   Confinement  therefore 
implies  that  v,"  «  v,,,  which  by  (33)  is  equivalent  to  requiring 
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23c.  ,  2 

B   :    2.8  webers/m 
z 


Figure   23 

Plasmoid  Flow  in  the  Curved  Solenoid, 
II.   Plasmoid  at  10  joules 
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2l^a. 
No  field 
(triple  shot) 


2[|b.  2 

B^:  1,14.  webers/in 

(double  shot) 


2l\.c, 
B^:  2.8  vjehera/m" 


Figure  21; 

Plasmoid  Flow  in  the  Curved  Solenoid. 
III.   Plasmoid  at  2.5  joules 
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that  V||  «  Ro)  =  (R/r)ra)  =  (R/r)yL  .   Since  v,,  =  3  v^  while 
R  »  r,  the  condition  is  clearly  satisfied:  had  the  plasmold  a 
sufficiently  low  density,  it  would  follow  the  field  lines. 
On  the  other  hand  a  sufficiently  dense  plasmold  will 
develop  a  polarization  electric  field,  as  the  result  of  charge 
separation  created  by  v ," ,  which  will  result  in  an  E  X  B  radial 
drift  of  almost  the  magnitude  needed  to  compensate  for  the 
field  curvature.  The  plasmold  will  thus  travel  in  a  nearly 
straight  line,  at  an  increasing  inclination  to  the  field.   A 
similar  effect  was  seen  in  the  experiments  of  Bostick,  where 
plasmold  velocity  was  perpendicular  to  the  field.   To  study  this 
(the  following  argument  was  kindly  suggested  by  Professor 
George  Schmidt  of  Stevens  Institute  of  Technology),  adopt  a 
rotating  coordinate  system  (traveling  with  the  plasmold)  of 
angular  velocity  v. /R.  The  charge  separation  drift,  v^" , 
creates  an  electric  displacement,  D,  which  Increases  at  the  rate 

D  =  qn7^  +  e^E,  (5^) 

with  q  the  charge  per  ion  and  n  the  nvmiber  density.   Now  assvune 
the  condition  (to  be  shortly  verified)  that 

€  »  €^  ,  (35) 

SO  that  the  second  term  to  the  right  in  {J>k)   becomes  negligible. 
Also,  since  the  polarizablllty,  e  -  e^,  of  the  plasmold  may  be 
taken  as""-^  Mn/B^  (with  M  the  ion  mass),  it  follows  from  (35) 
that 

€  =  Mn/B^.  (36) 
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Now,  the  radial  drift,  v  =  E  X  b/B^  will  increase  at  the  rate 


v^  =  E  X  b/B^  =  D  X  B/eB^  .  (37) 

If  one   combines   (33)   throxigh   (37),   he  obtains 

2 


v_   =   Vj"0O   =   v„    /R 


II, 
r  ~   'd 


which  implies  rectilinear  motion  in  the  laboratory  frame. 

However,  if  the  plasmoid  is  traveling  through  a  conducting 
background,  circulating  external  currents  are  created  (through 
the  background)  which  reduce  the  polarization  field  in  the 
laboratory  frame.   Indeed,  if  the  background  medium  remains  at 
rest,  leaving  the  plasmoid  moving  only  in  the  magnetic  field, 
the  plasmoid  will  now  follow  the  field  closely  (providing, 
further,  that  the  plasmoid  pressure  is  small  compared  to  B  /2m,  , 
the  p  «  1  condition). 

The  conditions,  then,  for  rectilinear  motion  are  a  vacuiwi 

region  surrounding  the  plasmoid  and  (35)  (which  may  be  seen 

2      P 
to  be  equivalent  to  Mnc  »  B  /2\i^).     For  the  present  experiment, 

taking  M-~6o  a.m.u.  and  B  =  2  weber/m^,  one  finds  that  (35)  and  > 

(36)  imply  that  n  »  3  x  10  ions/cm^  which  is  well-satisfied 

(Section  IIIB3).   However,  in  some  results,  e.g..  Figure  24,  an 

earlier  portion  of  the  plasmoid,  traveling  rectilinearly  into 

the  wall,  may  have  created  the  necessary  background  atmosphere 

for  permitting  later  portions  to  follow  the  field,  in  the  manner 

described  in  the  previous  paragraph.   The  results  of  Figure  22  and 

the  slight  effect  of  doubling  the  field  noted  in  Figure  23  probably 

relate  to  poor  satisfaction  of  the  3  condition. 
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APPENDIX  I 


Experimental  Equipment 


This  section  describes  the  overall  arrangement  of  experi- 
mental equipment.   The  Appendices  which  follow  will  provide 
details  on  the  various  components. 

Figure  25a.  shows  the  Differential  Solenoid  emplaced  upon 
the  plasma  source  capacitori  the  leads  to  the  search  coil  may  be 
seen.   Behind  is  the  solenoid  capacitor  bank  and  to  the  left  the 
ignitron  unit.   In  the  left  foreground  is  the  control  console 
Including  the  pulse  and  waveform  generators.  Vacu\am  equipment 
is  to  the  right,  off  the  picture.   Figure  25b.  is  a  closer  view. 
Here,  the  main  vacuum  system  is  in  the  background,  right.   In 
the  foreground  is  the  button  source  trigger  to  the  spark  gap  of 
the  main  (coaxial)  source,  together  with  the  vacuvmi  arm  (coming 
forward  and  up)  to  the  auxiliary  vacuum  system,  off  the  picture 
to  the  right . 

A  direct  view  of  the  trigger  system  is  given  in  Figure  26. 
The  cross  of  the  T  is  the  button  source  trigger,  with  capacitor 
at  right.   The  vacuum  line  comes  off  forward.   At  left  is  a 
top  view  of  the  capacitor  for  the  coaxial  plasma  source  which 
is  mounted  above  it. 

Figure  27  gives  two  views  of  the  Curved  Solenoid.   Here  the 
button  source  assembly  is  covered  with  photographer's  cloth  as 
in  an  actual  run  (to  block  spurious  light  to  the  camera).   The 
search  coil  is  clearly  seen. 
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The  entire  experimental  arrangement  Is  shown  schematically 
In  Figure  28.   The  main  power  supply,  control  unit,  Ignitron 
unit,  and  pulse  amplifier  were  designed  and  built  in  this  labora- 
tory and  are  described  in  Appendix  IV.  The  coaxial  plasma  source 
was  also  original  and  is  described  in  Appendix  II.   The  trigger 
("ion  gun"),  usually  a  button  or  rail  source,  is  described  in 
Appendix  III»   The  main  vacuum  system  (6")  was  a  commercially 
assembled  unit  (Kinney  Model  PW-6OO);  the  auxiliary  system  (2") 
was  laboratory  assembled  around  a  commercial  diffusion  pxamp 
(CVC  Model  MCF-60).   The  5  kV  supply  was  a  laboratory  assembled 
half-wave  rectifier  with  a  2  kVA  transformer,  and  some  filtering 
to  counteract  the  small  time  constant  associated  with  the  trigger 
capacitor.   The  12  kV  supply  was  a  commercial  30  kV,  2raA  unit 
(Beta,  Model  202/R).   The  main  capacitor  bank  was  described  in 
Section  IIC.   The  pulse  and  sawtooth  generators  were  commercial 
units  (Tektronix,  Models  I6I  and  I62,  respectively). 


-  88  - 


Figure  25 
Experimental  Arrangement  with  Differential  Solenoid 
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Figure   26 
Button  Source   Trigger 
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27a. 


27b. 


Figure  27 
Experimental  Arrangement  with  Curved  Solenoid 
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APPENDIX  II 
Coaxial  Plasma  Source 

The  main  source  for  the  experiments  was  a  pair  of 
coaxial,  brass  electrodes,  employing  the  "shaped  charge" 
geometry.   Design  details  are  seen  in  Figure  29.   The 
center  terminals  of  the  firing  capacitor  were  fastened 
through  a  cylindrical  brass  block,  out  of  which  projected 
the  bottom  electrode  of  the  trigger  spark  gap.   The  upper 
electrode  of  the  spark  gap  and  the  center  electrode  of 
the  plasma  source  were  connected  and  seated  in  Teflon-   The 
return  path  of  the  outer  plasma  electrode  was  along  a 
section  of  2"  brass  pipe  (appropriately  slotted)  which  was 
seated  on  a  circular  copper  plate  connected  to  the  ground 
terminals  of  the  capacitor. 
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APPENDIX  III 
Button  Source  Trigger 

The  trigger  for  the  vacuum  spark  gap  switch  of  Figure  29 
is  shown  In  Flg\ire  26,  Inserted  Into  the  cross-bar  of  the  T, 
with  Its  own  coaxial  capacitor  on  the  right.   A  detail  of  this 
arrangement  Is  shown  In  Figure  30  (where.  In  this  drawing,  the 
source  was  shown  coaxial).   The  capacitor  was  a  hollow  cylinder 
of  barlxmi  tltanate,  slivered  on  the  Inner  and  outer  surfaces. 
Trigger  was  air  gap  from  the  pulse  amplifier  (see  Appendix  IV). 
This  figure  shows  the  arrangement  which  was  actually  employed 
as  the  main  plasma  source  of  earlier  experiments;  when  the 
larger  coaxial  source  was  constructed  (Appendix  II),  a  button 
gun  was  substituted  on  this  device  and  it  was  employed  as  the 
spark  gap  trigger. 
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Figure  29 
Detail  of  Coaxial  Plasma  Source 
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APPENDIX  IV 
Laboratory  Designed  Circuits 

There  are  two  circuits  to  be  distinguished;  the  pulse 
amplifier,  and  the  main  power  supply  and  firing  system  for  the 
capacitor  bank  and  solenoid. 

A  rather  simple  pulse  amplifier  was  designed  to  magnify 
the  signal  from  the  pulse  generator  (see  Appendix  I)  to  an 
amplitude  sufficient  to  trigger  the  air  gap  of  the  barium 
titanate  capacitor  employed  with  the  plasma  trigger  (see 
Appendix  III).   The  circuit  is  shown  in  Figure  51.   It  was  found 
that  a  20  V  trigger  satisfactorily  tripped  the  2D21  at  about 
-12  V  bias.   The  resulting  +200  V  swing  was  adequate  to  trigger 
the  5C22. 

The  main  power  supply  for  charging  the  capacitor  bank  is 
shown  in  Figure  32.   The  plate  transformer  was  a  6  kVA  commercial 
item  (United  Transformer,  Model  CQ-309),  operated  on  two  phase 
(208  V)  with  one  leg  direct  and  one  fed  through  a  motor  driven 
2.8  kVA  variable  transformer  (Powerstat  Model  MB  I36).   Recti- 
fication was  provided  by  a  mercury  vapor  tube  bridge.   The 
variable  transformer  crank -up  time  was  5  sec.   Charging 
resistance  to  the  bank  (840  ijlP,  3  kV)  was  usually  2000  ohms 
(500  ohms  in  ignitron  unit). 

The  ignitron  unit  is  shown  in  Figure  33.   The  GL  5550  and 
the  GL  7171  have  maximum  pulse  current  ratings  of  900  amps  and 
30,000  amps,  respectively.   The  former  is  shown  here  triggering 
the  latter  at  3  kV  and  200  amps.   The  trigger  for  the  GL  5550 
was  taken  from  the  discharge  of  20  m.F  at  600  V,  as  shown,  the 
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capacitors  being  charged  by  a  variable  voltage  divider,  off  B  , 
that  could  be  adjusted  to  compensate  for  various  desired  firing 
voltages  on  the  main  capacitor  bank  (the  circuit  is  shown  here 
adjusted  for  3  kV).   The  0.001  \iF ,   10  \i\iF ,   and  0.1  tiF  capacitors 
were  employed  for  noise  suppression.   The  one  megohm  voltage 
divider  circuit  off  the  ignitor  was  used  occasionally  as  an 
oscilloscope  trigger.   The  label  "coil"  refers  to  the  main 
solenoid  to  be  energized.   The  lOK  resistance  was  employed  to 
protect  the  bank  from  isolation.  Buzzer  and  bell  were  warnings 
of  "ready"  and  "charging". 

The  control  unit  for  the  charging  and  firing  system  is 
shown  in  Figure  34.   This  unit  was  designed  to  perform  the 
following  functions;  provide  filament  and  control  power  and 
time  delay  for  the  filaments,  protect  against  the  variable  trans- 
former drive  motor  being  in  the  off -zero  position  or  lacking 
activating  power,  provide  a  signal  for  activating  the  main  power 
supply  to  begin  the  bank  charging  cycle,  provide  a  voltmeter 
sensing  unit  (meter  relay  circuit)  for  deactivating  the  power 
supply  and  firing  the  ignitron  unit  at  a  pre-set  voltage, 
provide  means  for  immediate  (or  emergency)  firing  of  the  ignitron 
unit  at  any  stage  of  the  charging  cycle  (if  the  voltage  is  too 
low,  the  ignitron  will  not  fire  and  the  bank  will  discharge 
through  the  lOK  resistor),  and  provide  primary  current  monitoring 
on  one  input  leg  of  the  power  supply.   There  are  also  provisions 
for  stopping  the  drive  motor  of  the  variable  transformer  at  any 
time  after  initiation  of  the  charge-fire  cycle,  and  for  extinguish- 
ing all  console  and  indicator  lights. 
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Details  of  the  meter  relay  circuit  are  shown  in  Figure 
35-   It  is  seen  that  the  VI -V2  circuit  is  not  closed  unless 
the  DC  circuit  for  powering  the  relay  is  operative  (i,e,, 
fail-safe).   The  reversing  switch  compensates  for  change  in 
bank  polarityo   The  "meter  hold"  switch  is  usually  left 
open  so  that  appropriate  relays  in  the  control  unit  can 
release  the  meter  (which  is  of  the  "locking"  variety)  after 
firing. 
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APPENDIX  V 
Photomultiplier  Experiment 

The  center-of-raass  and  thermal  velocity  results  were 
obtained  from  oscilloscope  traces  of  photomultiplier  (PM)  out- 
puts, the  PM  assemblies  being  stationed  at  three  positions  along 
the  experiment  tube,  as  described  in  Section  IIG. 

Each  assembly  consisted  of  a  6292  PM  (Tracerlab  Model  RLD-2 
Gamma  Scintillation  Spectrometer  Detector  with  the  sodivun  iodide 
crystal  removed),  powered  at  -2500  V  (Tracerlab  Model  RLS-l), 
A  three  foot  colliraating  tube  was  attached  to  the  front  of  the 
PM  with  a  thin  slit  in  the  cover  at  the  far  end,  which  was 

oriented  normal  to  the  experiment  axis.   Thus,  the  error  in 
position  at  each  PM  station  was  of  the  order  of  one  millimeter. 
The  sensitivity  of  the  6292  is  extremely  large j  consequently, 
several  neutral  density  filters  were  required  for  each  assembly, 
in  some  cases,  reducing  the  maximiom  light  intensity  by  four  to 
six  orders  of  magnitude  in  order  that  the  tube  not  be  overdriven. 
The  outputs  of  the  three  PM  assemblies  were  each  run  to  ground 
through  an  appropriate  dropping  resistor,  the  potentials  developed 
being  read  on  separate  oscilloscopes  (Tektronix,  Model  5^5,  with 
Type  53/5^K  Plug-In  Preamplifier).   Dropping  resistances  were 
in  the  range  of  500  to  1000  ohms,  with  current  maxima  of  0.3  to 
3  mA. 

Oscilloscope  triggering  was  arranged  in  the  following  fashion. 
Time  zero  was  taken  to  be  the  firing  of  the  trigger  spark  gap  on 
the  main  plasma  source.   A  929  phototube  was  focused  on  the  gapj 
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the  emitter  was  kept  at  -500  V  by  a  regulated  power  supply  while 
the  anode  was  taken  to  ground  through  an  appropriate  dropping 
resistance.   The  potential  developed  (negative)  triggered  the 
first  oscilloscope  (the  PM  assembly  closest  to  the  plasma 
source).   A  rather  curious  fact  emerged  in  connection  with 
designing  this  arrangement.   If  the  plasma  source  was  fired 
without  the  main  confining  field  being  activated  (the  beginning 
of  the  solenoid  was  a  foot  or  more  above  the  spark  gap)  the 
light  from  the  spark  gap  produced  a  triggering  signal  of  approxi- 
mately -10  V  across  350  ohms.   If  the  main  solenoid  was  activated 
also  (and  was  therefore  at  maximum  field  when  the  spark  gap 
was  energized)  the  spark  gap  light  was  considerably  diminished j 
a  dropping  resistance  of  3900  ohms  was  required  to  develop  a 
satisfactory  triggering  potential. 

The  second  oscilloscope  was  triggered  from  the  positive 
gate  of  the  first,  and  the  third  from  that  of  the  second.   It 
was  established  that  each  gate  Introduced  a  delay  of  approxi- 
mately  0.1  M-sec.   The  sweep  times  of  the  three  oscilloscopes 
had  been  synchronized  on  the  sweep  scales  actually  employed. 


-  106  - 


APPENDIX  VI 

Determination  of  Flux  Distribution  in 
Multiple  Circular  Coils  with  Axial  Ssntunetry 


In  Section  IIEl  it  was  indicated  that  a  circular  model 
with  axial  symmetry  could  be  employed  to  approximate  the  helical 
windings  (Differential  Solenoid)  employed  in  this  experiment. 
The  procedure  consists,  first,  of  approximating  each  turn  by  a 
circular  filament  and  reducing  the  model  to  two  dimensions o 
Then  the  field  created  by  the  filament  is  obtained  for  any 
position  on  the  plane,  and  thus  in  space,  as  a  function  of  the 
current  in  the  filament,  by  application  of  the  Biot-Savart  law 
(actually,  one  obtains  B  ,  B  ,  and  (^/2Tr,  a  stream  function). 
If  there  are  many  turns,  then  the  effect  of  all  the  filaments 
on  the  point  in  question  must  be  summed. 

Once  the  results  are  obtained  for  the  first  poimt  in  the 
(r  >  0,z)  half-plane,  one  can  then  move  in  the  z  (axial)  direc- 
tion a  small,  convenient  distance  and  then  proceed  to  search  in 
the  radial  direction  for  another  point  of  the  same  stream-function 
magnitude.   Having  found  the  point,  B  and  B  can  be  computed  again 
and  the  process  continued.   In  this  way,  by  choosing  various 
starting  positions  for  streamlines,  the  entire  field  can  be 
mapped  to  any  desired  degree  of  completeness. 

It  may  be  seen  that  this  procedure  has  wide  applicability. 
Consider  a  coil  with  varying  spacing  between  turns  and  varying 
radius,  or  a  group  of  colls  in  which,  in  addition,  the  current 
is  altered.   If  cylindrical  symmetry  obtains,  the  model  is 
applicable.   Even  more  complex  is  the  single  turn  of  heavy  copper, 
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perhaps  machined  in  such  a  way  as  to  alter  the  current  distribu- 
tion, which  is  particularly  popular  in  some  "mirror  machine" 
designs.   By  dividing  the  cross-section  of  the  turn  into 
Individual  filaments  with  appropriate  currents,  the  model 
results o   The  accuracy  of  the  model  is  improved  by  increasing 
the  number  of  subdivisions. 

This  problem  is  ideal  for  digital  machine  computation  and 
was  Fortran-programmed  for  the  IBM  704,   Input  consisted  of  th« 
coordinates  of  the  intersection  of  each  filament  with  the  r-z 
half -plane,  together  with  the  current,  and  the  coordinates  of 
one  or  more  starting  positions  of  streamlines.   One  specified 
also  the  length  of  the  Jump  in  both  the  z  and  r  directions  for 
searching  the  next  point,  the  degree  of  permitted  error  in  the 
magnitude  of  the  stream  function,  and  a  maximum  r  and  z  as 
cut-off  for  each  streamline.   There  was  provision  for  relocatiiig 
a  given  streamline  further  down  the  axis  which  may  have  become 
cut  off  in  r  as  a  result  of  a  gap  in  the  coil  winding. 

For  those  who  may  desire  to  employ  this  code,  located  at 
the  Institute,  it  is  identified  as  NYU  22.?.   Input  symbols  ares 
C,  current  in  a  current  filament i  PSI,  the  stream  function 
magnitude,  equal  to  ^  /2ir|  RHO,  the  radial  location  of  the 
intersection  of  a  current  filament  in  the  half -plane i  ZETA,  the 
axial  location  of  a  current  filament i  RS,  the  starting  radial 
position  for  the  plot  of  a  streamline  at  the  first  axial  value 
of  ZSj  ZS,  the  starting  axial  position  for  the  plot  of  a  stream- 
line and  any  subsequent  axial  position  at  which  it  is  to  be 
relocated  I  NCOIL,  the  nxamber  of  current  filaments  being  employed 
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in  the  model  I  NLINE,  the  niimber  of  streamlines  to  be  plotted  i 
NPOSN,  the  nximber  of  additional  axial  positions  (ZS),  plus  one, 
for  which  streamlines  already  requested  are  to  be  relocated 
(assuming  that  some  streamlines  have  been  cut  off  by  having 
acquired  too  large  a  radial  component j  if  a  streamline  was  not 
cut  off  then  the  machine  will  ignore  the  order  to  recompute  it  at 
the  new  ZS)j  RT,  the  radial  dimension  at  each  value  of  ZS  beyond 
which  the  machine  will  not  attempt  to  relocate  a  streamline i  H, 
the  trial  step  in  the  axial  direction  for  locating  the  next  point 
on  the  same  streamline j  DR,  the  trial  step  in  the  radial  direction 
(after  having  moved  H)  for  locating  a  point  the  streamline  magni- 
tude of  which  will  differ  from  the  previous  point  an  amount  less  than 
ERRORi  ZMIN,  ZMAX,  line-plot  cut-offs  in  the  axial  direction; 
RMAX,  same  for  the  radial  direction. 

In  the  printed  output  each  page  (or  group  of  pages)  will 
contain  the  data  for  one  streamline,  headed  by  its  magnitude, 
and  including  every  point  on  the  line,  together  with  the  values, 
for  that  point,  of  the  magnetic  field  intensity  (B)  and  its  radial 
and  axial  components  (BR,  BZ).  The  numerical  answers  are  given  in 
field  intensity,  not  induction  (despite  the  use  of  "B" ) ,  and 
if  linear  dimensions  are  given  in  meters  and  current  in  amperes, 
the  results  will  be  unrationalized  MKS  (i.e.,  \i     =  10"  ). 
Similarly,  $  =  2TruQPSI. 

There  are  some  words  of  caution  appropriate  to  this  machine 
code.   As  was  previously  noted,  the  model  becomes  decrea singly 
applicable  in  regions  close  to  the  current  sheath  (i.e.,  the  coils') 
if  the  actual  winding  Involves  helical  rather  than  cylindrical 
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symmetry.   Also,  In  obtaining  line  plots  close  to  colls,  small 
test  jumps  are  recommended,  of  magnitude  significantly  smaller 
than  the  distance  to  the  nearest  coll  (filament),  else  the 
machine  will  become  "lost"  and  not  find  the  next  point.   Since 
the  mathematical  expressions  employed  Involve  elliptic  integrals, 
care  must  be  taken  in  this  region  (near  the  colls)  that  the 
machine  sub-routine  for  the  integral  is  employing  a  sufficient 
number  of  terms  (or  that  there  is  some  alternate  provision), 
as  a  modulus  close  to  unity  (which  is  the  result  of  a  close 
approach  to  a  filament)  produces  an  approach  to  infinity  in 
Integrals  of  the  first  kind.   Finally,  note  that  the  value  of 
B  on  the  axis  may  not  be  obtained  since  the  stream  function  is 
zero  there  and,  thus,  smaller  than  the  smallest  permissible 
error  (in  fact,  the  program  will  skip  such  a  request). 
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